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Abstract. Recent months have witnessed dramatic progress in our understanding of 
' short 7-ray burst (SGRB) sources. There is now general agreement that SGRBs - or 

. at least a substantial subset of them - are capable of producing directed outflows of 

I relativistic matter with a kinetic luminosity exceeding by many millions that of active 

OO ' galactic nuclei. Given the twin requirements of energy and compactness, it is widely 

. believed that SGRB activity is ultimately ascribable to a modest fraction of a solar 

J — I mass of gas accreting onto a stellar mass black hole or to a precursor stage whose 

■ inevitable end point is a stellar mass black hole. Astrophysical scenarios involving the 
fH , violent birth of a rapidly rotating neutron star, or an accreting black hole in a merging 
Qh' compact binary driven by gravitational wave emission are reviewed, along with other 
Q . possible alternatives (collisions or collapse of compact objects). If a black hole lies 
| H I at the center of this activity, then the fundamental pathways through which mass, 
^ ' angular momentum and energy can flow around and away from it play a key role in 
. . , understanding how these prime movers can form coUimated, relativistic outflows. Flow 

. ^ ' patterns near black holes accreting matter in the hypercritical regime, where photons 

^ ■ are unable to provide cooling, but neutrinos do so efficiently, are discussed in detail, 

^ I and we believe that they offer the best hope of understanding the central engine. 

- - ' On the other hand, statistical investigations of SGRB niches also furnish valuable 

information on their nature and evolutionary behavior. The formation of particular 
kinds of progenitor sources appears to be correlated with environmental effects and 
cosmic epoch. In addition, there is now compelling evidence for the continuous fueling 
of SGRB sources. We suggest here that the observed late flaring activity could be 
due to a secondary accretion episode induced by the delayed fall back of material 
dynamically stripped from a compact object during a merger or collision. Some 
important unresolved questions are identified, along with the types of observation 
that would discriminate among the various models. Many of the observed properties 
can be understood as resulting from outflows driven by hyperaccreting black holes 
and subsequently collimated into a pair of anti-parallel jets. It is likely that most 
of the radiation we receive is reprocessed by matter quite distant to the black hole; 
SGRB jets, if powered by the hole itself, may therefore be one of the few observable 
consequences of how flows near nuclear density behave under the influence of strong 
gravitational fields. 
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1. Introduction 

1.1. Prologue 

In the sections which follow, we shall be concerned predominantly with the theory 
of short 7-ray burst^. If the concepts there proposed are indeed relevant to an 
understanding of the nature of these sources, then their existence becomes inextricably 
linked to the metabolic pathways through which gravity, spin, and energy can combine 
to form coUimated, ultra relativistic outflows. These threads are few and fragile, as we 
are still wrestling with understanding non-relativistic processes, most notably those 
associated with the electromagnetic field and gas dynamics. If we are to improve 
our picture-making we must make more and stronger ties to physical theory. But in 
reconstructing the creature, we must be guided by our eyes and their extensions. In 
this introductory chapter we have therefore attempted to briefly summarize the observed 
properties of these ultra-energetic phenomena|§|. There are five sections: §1.2| gives a brief 
account of their history from birth to present-age; §1.31 is devoted to their metabolism 
- in other words, to their gross energetics, spectra and time variability; §1.41 describes 
the attributes of the afterglow signals, which, as fading beacons, mark the location 
of the fiery and brief 7-ray event. These afterglows in turn enable the measurement 
of redshift distances, the identification of host galaxies at cosmological distances, and 
provide evidence that many short 7-ray bursts are associated with old stellar populations 
and possibly with no bright supernova. These threads will be woven in §1.51 Finally, 
§1.61 gives a compendium of the observational facts. 



1.2. Burst of Progress 

The manifestations of SGRB activity are extremely diverse. SGRBs are observed 
throughout the whole electromagnetic spectrum, from GHz radio waves to 10 MeV 
7-rays, but until recently, they were known predominantly as bursts of 7-rays, largely 
devoid of any observable traces at any other wavelengths. 

Before 2005, most of what we knew about SGRBs was based on observations 
from the Burst and Transient Source Experiment (BATSE) on board the Compton 
Gamma Ray Observatory, whose results have been summarized by Fishman & Meegan 
[2]. BATSE, which measured about 3000 events, detected approximately one burst 
on a typical day. While they are on, they outshine every other source in the 7-ray 
sky, including the sun. Although each is unique, the bursts fall into one of two rough 
categories. Bursts that last less than two seconds are classified as short, and those that 
last longer - the majority - as long [3]. The two categories differ spectroscopically, with 

I The literature on this subject has become quite large, and to keep the references to a manageable 
size, we have in general referred to the most recent comprehensive article in a given topic. In particular, 
we give only representative references for those topics on which helpful reviews articles already exist. 
It is thus possible for the reader to trace the early work on which some of these conclusions are based. 
We apologize to those colleagues whose work is either omitted or not fully represented. 
§ The reader is referred to for an an excellent review of the observations. 
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Figure 1. Hardness ratio (i.e., ratio of total counts in the 100-300 keV range, to 
those in the 50-100 keV range) versus duration for BATSE bursts. The hardness 
ratio is a measure of the shape of the spectrum: larger values correspond to harder 
spectra. Different line levels denote number density contours. The bimodality of the 
distribution of their duration is confirmed by the associated spectral shape. 

short bursts having relatively more high-energy 7-rays than long bursts do. Figure [1] 
shows the hardness ratio as a function of the duration of the emission. It is a measure 
of the slope of the spectrum, where larger values mean that the flux at high energies 
dominates. 

Arguably the most important result from BATSE concerned the spatial distribution 
of bursts. Both long and short events occur isotropically - that is, evenly over the entire 
sky with no dipole and quadrupole components, suggesting a cosmological distribution. 
This finding cast doubt on the prevailing wisdom, which held that bursts came from 
sources within the Milky Way. The uniform distribution instead led most astronomers 
to conclude that the instruments were picking up some kind of cosmological event. 
Unfortunately, 7-rays alone did not provide enough information to settle the question 
for sure. The detection of radiation from bursts at other wavelengths would turn out 
to be essential. Visible light, for example, could reveal the galaxies in which the bursts 
took place, allowing their distances to be measured. Attempts were made to detect 
these burst counterparts, but they proved fruitless. 

Observations of burst counterparts ^ were restricted to the class of long duration 
bursts [Jj] until, in 2005, the Swift spacecraft succeeded in obtaining high-resolution X-ray 
images [5] [6] of the fading afterglow of GRB 050509B - so named because it occurred 
on May 9, 2005. This detection, followed by a number of others at an approximate rate 
of 10 per year, led to accurate positions, which allowed the detection and follow-up of 
the afterglows at optical and longer wavelengths [3, [HI [9] . This paved the way for the 
measurement of redshift distances, the identification of candidate host galaxies, and the 
confirmation that they were at cosmological distances P [TOl [TT|, |5l El [121 [13] ■ Swift is 

11 because BeppoSAX is mainly sensitive to bursts longer than about 5 to 10 s. 
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Figure 2. Representative spectra lyF^ (x v^N{v) of various SGRBs. Tlie SGRB 
spectrum is non thermal, the number of photons varying typically as N{v) oc , 
where a ~ 1 at low energies changes to a ~ 2 to 3 above a photon energy ^ 1 MeV. 



equipped with 7-ray, X-ray and optical detectors for on-board follow-up, and capable of 
relaying to the ground arc-second quality burst coordinates within less than a minute 
from the burst trigger, allowing even mid-size ground-based telescopes to obtain prompt 
spectra and redshifts. 



1.3. Metabolics 

SGRBs are brief flashes of radiation at soft and hard 7-ray energies that display a wide 
variety of time histories. They were first detected at soft 7-ray energies with wide field- 
of-view instruments, with peak soft 7-ray fluxes reaching hundreds of photons cm~^ 
s~^ in rare cases. The BATSE instrument was sensitive in the 50-300 keV band, and 
provided the most extensive data base of SGRB observations during the prompt phase. 

SGRBs typically show a very hard spectrum in the soft to hard 7-ray regime. The 
photon index breaks from ^ —1 at energies Ep^ < 100 keV, to a —2 to —3 spectrum 
at Epii > several hundred keV [H]. Consequently, the peak photon energies, i?pk, of 
the time-averaged uF^y spectra of BATSE SGRBs are typically found in the 500 keV 
- several MeV range [151 ESI [U] • The general trend is that the spectrum softens, and 
Ep]^ decreases, with time. More precise statements must, however, wait for larger area 
detectors. 

In Figure [2] representative spectra are plotted, in the conventional coordinates 
u and z/Fjy, the energy radiated per logarithmic frequency interval. Some obvious 
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Figure 3. BATSE lightcurves of various individual bursts. The y axis is the photon 
count rate in the 0.05 to 0.5 MeV range; the x axis is the time in seconds since the 
burst trigger. Both before and after the burst trigger, no 7-rays are detected, above 
background, from the same direction. 



points should be emphasized. We measure directly only the specific luminosity 
D^vljj = {1/4:'k)vLi, (the energy radiated in the direction of the earth per second 
per steradian per logarithmic frequency interval by a source at luminosity distance 
D), and its dimensionless distance-independent ratios between two frequencies (BATSE 
triggers, for example, are based on the count rate between 50 keV and SOOkeV). The 
apparent bolometric luminosity 47r D^vlyd{hiv) may be quite different from the true 
bolometric luminosity /q°° D^I^dudfl if the source is not isotropic. The GRB spectra 
shown in Figure [2] are those of GRB 910609, GRB 970704, and GRB 931101 which were 
observed by BATSE. The time integrated spectrum on those detectors ranges from 25 
keV to lOMeV [TT]. 

A typical SGRB - if there is such a thing - lasts for a fraction of a second. Observed 
durations vary, however, by three orders of magnitude, from several milliseconds [2] to 
a few seconds [3]. The shortest BATSE burst had a duration of 5ms with a 0.2ms 
structure [18]. Similarly to long bursts, SGRBs have complicated and irregular time 
profiles which vary drastically from one burst to another [19]. They range from smooth, 
fast rise and quasi-exponential decay, through curves with several peaks, to variable 
curves with many peaks. Various profiles, selected from the BATSE catalog, are shown 
in Figure [3l 

The duration of a SGRB is defined by the time during which the middle 50% 
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Figure 4. Apparent isotropic 7-ray energy as a function of redshift. The energy is 
calculated, assuming isotropic emission, for 6 SGRBs with estimated redshifts. For 
comparison, the isotropic 7-ray energies for 22 long GRBs are also plotted using the 
compilation of [2T] . 



(tso) or 90% (tgo) of the counts above background are measured. A bimodal duration 
distribution is measured, irrespective of whether the ^50 or tgQ durations are considered 
[3]. About two-thirds of BATSE GRBs are long-duration GRBs with tgo > 2 s, with the 
remainder comprising the SGRBs. The integral size distribution of BATSE SGRBs in 
terms of peak flux 0p is very flat below ~ 1 ph cm~^ s~^, and becomes steeper than the 
expected power law with index —3/2 of a Euclidean distribution of sources, at 0p > 5 ph 
cm~^ s~^ (see Figure 12 in [2]). This follows from a cosmological origin of GRB sources, 
with the decline in the number of faint bursts due to cosmic expansion. Follow-up X-ray 
observations with Swift and HETE-IIhave permitted redshift determinations that firmly 
establish the distance scale to the sources of SGRBs [SniElElElIIIlIISlliaESlEniEllEH]. 

The redshifts of about half a dozen SGRBs are now known (late 2006), with the 
median {z) ~ 0.3 and the largest spectroscopically mieiied redshift at z = 0.71. It should 
be noted that the inference of redshift distances relies at present on the statistical 
connection to a putative host galaxy and spectroscopy of the host, and that none 
are based on the absorption-line systems seen in the spectra of the afterglows. The 
corresponding distances imply apparent isotropic 7-ray energy releases -E-y^iso ~ 10^^- 
10^^ {Qj / Atc) ergs, where is the solid angle into which the 7-rays are beamed 
(Figure Hj). For a solar-mass object, this implies that an unusually large fraction of the 
energy is converted into 7-ray photon energy. This spread in the inferred luminosities 
obtained under the assumption of isotropic emission may be reduced if most GRB 
outflows are jet-like. A beamed jet would alleviate the energy requirements, and some 
observational evidence does suggest the presence of a jet (see §1.4p . 
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Figure 5. Swift lightcurves of various individual bursts: GRB 050724, GRB 051210, 
and GRB 051221A. The y axis is the flux in the 2 to 10 keV range; the x axis is the 
time in seconds since the burst trigger. 

1.4- A Warm Afterglow 

Among the first SGRBs localized by Swift was GRB 050509B. The satellite slewed 
rapidly in the direction of the burst, enabling its narrow- field X-ray telescope to 
pinpoint to within 10 arc seconds a source of faint X-ray emission less than a minute 
after the trigger. The source decayed with a power-law behaviour, (px oc t^, with 
X ~ —1.3 [5l[6]. Despite intense follow-up in the optical and radio, no counterparts were 
discovered [H [2Sl E] • However, the lack of a detectable afterglow at other wavelengths 
is not surprising considering, at face value, the existing GRB afterglow theory [27]. It 
took about two months before the observation of a second short burst, GRB 050709, 
pinpointed by HETE-II [13]. Other detections have followed, at an approximate rate of 
10 per year, and permitted the observation and follow-up of afterglows at optical and 
longer wavelengths. GRB 050709 was the first SGRB from which an optical counterpart 
was observed [71 E], and GRB 050724 was the first SGRB for which a radio afterglow 
was measured [9]. GRB 051221A is probably the best-sampled SGRB to date [22] . 

Of the few bursts localized by the Swift XRT, only four were bright enough to 
permit detailed study: GRB 050724, GRB 051221A, GRB 051210, and GRB 051227 
[281 [291 [30] . Various X-ray light curves selected from the Swift catalog, are shown in 
Figure [51 Chandra follow-up observations obtained for GRB 050724 and GRB 051221A 
provide the most constraining observations for SGRB jets to date [28l[29j. The lack of an 
observed early downturn in the X-ray light curves has been interpreted as supporting 
a mild degree of anisotropy in SGRBs [29l [28l [22l [31|, [32] • Anisotropy in the burst 
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outflow and emission affects the light curve at the time when the inverse of the bulk 
Lorentz factor equals the opening angle of the outflow. If the critical Lorentz factor 
is less than 3 or so (i.e. the opening angle exceeds 20°) such a transition might be 
masked by the change from ultra relativistic to mildly relativistic flow. It would then 
be generically difficult to limit the late-time afterglow opening angle in this way if it 
exceeds 20°. Since some afterglow light curves are unbroken power laws for over 20 days 
(e.g., GRB 050724), if the energy input were indeed just a simple impulsive shell the 
opening angle of the late-time afterglow at long wavelengths would probably be greater 
than 25°, i.e. ^lopt > 0.1 [28]. However, even this still means that the energy estimates 
from the afterglow assuming isotropy could be > 10 times too high. The beaming angle 
for the 7-ray emission, could be smaller, and is much harder to constrain directly. 

An observation that attracted much attention was the discovery of long duration 
(~ 100 s) X-ray flares following several SGRBs [131 E] after a delay of ~ 30 s (see, e.g., 
GRB 050724 in Figure [5]). There is also independent support that X-ray emission on 
these time scales is detected when lightcurves of many bursts are stacked [331 El] • These 
observations may indicate that some sources display continued activity (at a variable 
level) over a period of minutes [HI [32]. Additional structure in the light curves has 
also emerged from a continued analysis of some of these objects down to the faintest 
X-ray flux levels. In some bursts (e.g., GRBs 050724 and 051221A), the X-ray light 
curve exhibits evidence for a late "energy refreshment" to the blast wave, on time scales 
comparable to the afterglow time scale. A reason for the late rise would be, say, if slower 
moving ejecta catches up with the main shock front injecting a substantial amount of 
energy and momentum [351 [36]. However, there are other mechanisms that can produce 
rising afterglow fluxes [37l [38] . 

1.5. Galactic Hosts, Supernova Family Ties and Cosmological Setting 

1.5.1. Demography Starting with the detection of GRB 050509B, a growing body 
of evidence has suggested that SGRBs are associated with an older and lower-redshift 
galactic population than long GRBs and, in a few cases, with large (> 10 kpc) projected 
offsets from the centers of their putative host galaxies (Table [H [39j). 



Table 1. 


Basic properties of the host galaxies of SGRBs (adapted from 


[39]). 


GRB 


z 


host 


in a galaxy 


SFR 


Metallicity 






classification 


cluster 


(Mq yr-i) 




050509B 


0.225 


E 


Yes 


<0.1 


-1 


050709 


0.160 


Irr/late-type dwarf 


No 


> 0.3 


0.25 


050724 


0.258 


early (E+SO) 


No 


< 0.05 


0.2 


050813 


0.7227(1.8) 


E? 


Likely 


<0.2 


-1 


051221A 


0.5459 


late-type dwarf? 


7 


~1.5 


-1 


060502B 


0.287 


E? 


? 


<0.4 


-0.2 



The discovery of GRB 050509B and its fading X-ray afterglow [3] led to the first 
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redshift and host galaxy association [6J for a SGRB. A chance association with such 
a galaxy was deemed unlikely even under conservative assumptions and stood in stark 
contrast with the lines-of-sight of long GRBs, for which no association with an early-type 
host was ever made. Further Swift and HE TE- II detections of SGRBs have continued to 
support this hypothesis, though SGRBs are not ubiquitously found at large offsets and 
associated with early-type galaxies. GRB 050724 [El [231 HQ] , GRB 050813 [12], and 
GRB 060502B like 050509B, were found to be in close association with old, red galaxies. 
GRB 050724 had optical and radio afterglow emission that pinpointed its location to 
be within its red host, making the association completely unambiguous, though the 
association of GRB 050813 and GRB 060502B with any single host remains tentative 
|41j . The absence of observable Ha and [O II] emission constrains the unobscured star 
formation rates in these galaxies to < O.2M0 yr~^, and the lack of Balmer absorption 
lines implies that the last significant star forming event occurred > 1 Gyr ago. Based on 
positions of the afterglows, two bursts (050509B and 050813) are very likely associated 
with clusters of galaxies [l2j . 

Not all hosts lack active star formation: GRB 050709 [El [3 [TOl [8] and GRB 
05 1221 A [22] both had optical afterglows and were associated with galaxies showing 
evidence for current, albeit low, star formation. The host of GRB 051221A, moreover, 
exhibits evidence for an evolved stellar population. Despite the availability of both 
X-ray and optical afterglow locations, no nearby host has successfully been identified 
for either GRB 060121 or GRB 060313. These observations indicate that these SGRBs 
occurred during the past ~ 7 Gyr of the universe {z < 1) in galaxies with diverse 
physical characteristics. 

In contrast to what is found for SGRBs, all of the confirmed long GRB host galaxies 
are actively forming stars with integrated, unobscured SFRs ~ 1— IOMq yr~^ [ ^ H5 | H6]. 
These host galaxies have small stellar masses and bluer colors than present-day spiral 
galaxies (suggesting a low metallicity; [E]). The ages implied for the long-soft GRBs 
are estimated to be < 0.2 Gyr \^^, which is significantly younger than the minimum 
ages derived for the early-type galaxies found to be associated with SGRBs. The cluster 
environments of at least two SGRBs contrast strikingly with the observation that no 
well- localized long-soft GRB has yet been associated with a cluster [48] . 

On the whole, the hosts of SGRBs, and by extension the progenitors, are not 
drawn from the same parent population of long GRBs. SGRBs appear to be more 
diffusely positioned around galaxies, and their associated hosts contain a generally older 
population of stars. 

1.5.2. Soft Gamma-Ray Repeaters in Nearby Galaxies It has been noted [l9l [50] that 
the giant flare (GF) observed from the putative galactic magnetar source SGR1806-20 
in December 2004 [l9l [50] could have looked like a classical SGRB had it occurred much 
farther away, thus making the tell-tale periodic signal characteristic of the neutron star 
rotation in the fading emission undetectable. The two previously recorded GFs of this 
type, one each from SGR 0520-66 on 5 March 1979 [51] and SGR 1900 + 14 on 27 August 
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1998 [52j, would have been detectable by existing instruments only out to ~ 8 Mpc, 
and it was therefore not previously thought that they could be the source of SGRBs. 
The main spike of the 27 December event would have resembled a short, hard GRB if 
it had occurred within ~ 40 Mpc, a distance scale encompassing the Virgo cluster |50j . 
However, the paucity of observed giant flares in our own Galaxy has so far precluded 
observationally based determinations of either their luminosity function or their rate. 

In the magnetar model, SGRs are isolated neutron stars with teragauss exterior 
magnetic fields and even stronger fields within [53], making them the most magnetized 
objects in the Universe. The large-scale reorganization of the magnetic field is thought 
to produce the observed GFs. The formation rate of magnetars is expected to track that 
of stars, which is ~ 0.013 Mq Mpc"^ yr-i in our Galaxy [50]. This suggests that BATSE 
would have triggered on such events as SGRBs at a rate of 30(Ai'Gai/0.01 yr~^) yr~^, 
compared with the all-sky BATSE rate of about 150 yr~^. Here A^^cai is the average 
rate of GFs in the Galaxy similar to the 27 December event. The observed isotropic 
distribution of short BATSE GRBs on the sky and the lack of excess events from the 
direction of the Virgo cluster suggests that only a small fraction, < 5%, of these events 
can be SGR GFs within 40 Mpc, implying that A^^g < 3 x 10^'^ yr^^ on average for a 
Galaxy like our own [50] . 

Before Swift detected SGRBs, searches for nearby galaxies within narrow Inter 
Planetary network (IPN) error boxes revealed already that only up to ~ 15% of them 
could be accounted for by magnetars capable of producing GFs [54j. Finally, a search 
for SGRBs with spectral characteristics similar to that of GFs (thermal spectra with 
kT ~ 100 keV) in the BATSE catalog has concluded that a small fraction (up to a few 
percent) of them could have originated from magnetars in nearby galaxies |55j . 

Magnetars are thought to be formed during core-collapse events inside massive 
stars, and because of their relatively short lifetimes as observable sources, ~ 10^ yr, 
would naturally be located in predominantly star-forming galaxies, while essentially 
none should be seen in ellipticals. One possible distinction of these from the classic 
SGRB population may well come from radio observations, because their radio afterglows 
should not be detectable beyond ~ 1 Mpc [56l |57]. The fraction of SGR events among 
what are now classified as short GRBs may not be dominant, but it should be detectable 
and can be tested with future Swift observations. 

1.5.3. Supernova Partnership Current observational limits [l3l[6l[26] indicate that any 
supernova-like event accompanying SGRBs would have to be over 50 times fainter than 
normal Type la SNe or Type Ic hypernovae, 5 times fainter than the faintest known 
la or Ic SNe, and fainter than the faintest known Type II SNe. These limits strongly 
constrain progenitor models for SGRBs. 

The hmits derived for GRB 050709 (filled triangles) and GRB 050509B (empty 
triangles) are plotted in Figure [6] along with two SN hght curves as they would appear 
at z = 0.225. The Type Ic SNe plotted are the very energetic Type Ic SN 1998bw 
associated with the long GRB 980425 [58] and the faint, fast-rise Type Ic SN 19941 j59]. 
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Figure 6. Magnitudes and upper limits of SGRB afterglows, GRB 050709 (filled 
triangles) and GRB 050509B (empty triangles). For comparison, two Type Ic 
supernovae are also shown: SN 1998bw was associated with the long GRB 980425. 



The absence of a SN severely rules out models predicting a normal SN la associated 
with SGRBs. Likewise, current observations disfavor a progenitor uniquely associated 
with massive stellar collapse. Observations of long GRBs at 2; < 0.7 are consistent with 
all having SN bumps [60]. The situation is less clear at z < 0.4 where all but two long 
GRBs have had SN features: GRB 980425 p]; GRB 031203 [6T]; GRB 030329 [62j; 
GRB 011121 [63J; but no SN was detected in GRBs 060505 and 060614 down to a flux 
limit at least hundreds of times fainter than SN 1998bw [64] . 

1.5.4- Lifetimes The distribution of time delays between progenitor formation and 
explosion is not yet well understood. It can be constrained using, for example, the 
SGRB rate as a function of redshift [651 ES, ETJ [39]. If SGRB explosions lag star 
formation by some considerable amount of time, then the intrinsic redshift distribution 
of SGRBs would be significantly skewed to lower z when compared to the universal star 
formation rate (SFR). A cursory comparison of the redshift distribution of SGRBs with 
the SFR reveals what appears to be a significant time delay of a few Gyr. It should be 
noted that there are inherent biases in the discovery of a GRB at a given redshift which 
are often difficult to quantify, such as complex trigger efficiencies and non-detections. 

An alternative and perhaps less restrictive approach may be to use the rates of 
SGRBs in different types of galaxies [681 [69]. On average, early-type galaxies have their 
stars formed earlier than late-type galaxies, and this difference, together with the time 
delay between progenitor formation and SGRB outburst, inevitably leads to different 
burst rates in the two types of galaxies. For instance, the morphological types for the 
bursts in Tabled] reflect a higher incidence of early-type galaxies than type la supernovae 
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and suggest progenitor lifetimes significantly exceeding 6 Gyr [51] • One difficulty with 
this idea is that star formation processes in these two types of galaxies might not be 
identical. For example, elliptical galaxies can form by the merging of two gas-rich 
galaxies. Many globular clusters can form in the merging process, which could enhance, 
for example, the fraction of binary progenitors and also change the lifetime distribution. 

Obviously, these estimates are only sketchy and should be taken as an order of 
magnitude estimate at present. However, it should improve as more hosts are detected. 
On the other hand, a large progenitor lifetime would help explain the apparent high in- 
cidence of galaxy cluster membership. This can be naturally explained by the fact that 
galaxies in over dense regions form earlier in hierarchical cosmologies and thus make a 
substantial contribution to the local stellar mass inventory |69]. Detailed observations 
of the astrophysics of individual GRB host galaxies may thus be essential before strin- 
gent constraints on the lifetime of SGRB progenitors can be placed. If confirmed with 
further host observations, this tendency of SGRB progenitors to be relatively old can 
help differentiate between various ways of forming a SGRB. 

This concludes our compendium of the facts. For ease of reference in the chapters 
that follow, they have been assembled here with a minimum of speculative interpretation. 

1.6. Setting the Stage 

SGRBs have been observed assiduously throughout the electromagnetic spectrum 
only recently and although we now know much about their collective and individual 
properties, we are still long way from being sure how they operate. 

These SGRB sources involve energies that can exceed 10^° ergs, the mass equivalent 
of 1/10,000 of a sun. Compared with the size of the sun, the seat of this activity is 
extraordinarily compact, as indicated by rapid variability of the radiation flux on time 
scales as short as milliseconds. It is unlikely that mass can be converted into energy 
with better than a few (up to ten) percent efficiency; therefore, the more powerful SGRB 
sources must "process" upwards of 10~^Mq through a region which is not much larger 
than the size of a neutron star (NS) or a stellar mass black hole (BH). No other entity 
can convert mass to energy with such a high efficiency, or within such a small volume. 

Well-known arguments connected with opacity, variability time scales and so forth 
require highly relativistic and variable outflow [TOl E2l EH] . Best-guess numbers are 
Lorentz factors F in the range 10^ to 10'^, allowing rapidly-variable emission to occur 
at radii in the range 10^^ to 10^^ cm [TH [751 EHl [ZZj- Because the emitting region must 
be several powers of ten larger than the compact object that acts as trigger, there is 
a further physical requirement: the original energy outflowing in a wind would, after 
expansion, be transformed into bulk kinetic energy. This energy cannot be efficiently 
radiated as 7-rays unless it is re-randomized. The emitted energy is an observable 
diagnostic of the microphysical processes of particle acceleration and cooling occurring 
within the bulk flow [78|, [79l |80]. It is beyond the scope of this article to describe 
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the properties of these outflows and the physical processes occurring within, and we 
shall confine our attention to properties directly relevant to their origins (An excellent 
account of the structure and energetics of SGRB outflows is given in [Ij.). This review 
is therefore not complete and in its emphasis reflects the biases of the authors. 

Based on current models, the simplest hypothesis - that the afterglow is due to a 
relativistic expanding blast wave - seems to agree with the present data [271 [6l [22]. 
The complex time-structure of some bursts suggests that the central engine may remain 
active for up to 100 seconds [IHl IHll E]. However, at much later times all memory of 
the initial time-structure would be lost: essentially all that matters is how much energy 
and momentum has been injected, its distribution in angle and velocity. However we 
can at present only infer the energy per solid angle; there are reasons to suspect that 
the afterglow is not too narrowly beamed; on the other hand the constraints on the 
angle-integrated 7-ray energy are not strong [291 EHl [22], [3ll [32] . 

As regards the trigger, there remain a number of key questions. What is the iden- 
tity of their progenitors? What is the nature of the triggering mechanism, the transport 
of the energy and the time scales involved? Does it involve a black hole orbited by a 
dense torus? And, if so, can we decide between the various alternative ways of forming 
it? The presence of SGRBs in old stellar populations helps rule out a source uniquely 
associated with recent star formation, while the lack of an accompanying supernova is 
strong evidence against a core-collapse origin. There is now a stronger motivation to 
develop models in fuller detail. This article outlines some of these issues. 
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2. Basic Ingredients 

In this section, we present a partial summary of some general ways in which gravity, 
angular momentum and the electromagnetic field can couple to power ultra relativistic 
outfiows, along with a review of the most popular current models for the central source. 
There are four sections: §2.11 gives a brief account of the arguments in favor of a 
gravitationally fueled origin; §2.21 describes the attributes of the most widely favored 
and conventional progenitors; the various modes of energy extraction from such systems 
are then discussed in §2.31 finally, §2.41 gives a compendium of the types of observation 
that might help discriminate among the various progenitor models. 

2.1. General Gonsiderations 

Shortly after the discovery of quasars in 1963, it was suggested that accretion of gas onto 
a compact massive body was responsible for their enormous energy output [821 [83], which 
can be hundreds of times larger than that of entire galaxies. Such a body is essentially 
a very efficient converter of gravitational binding energy into radiation. The deeper the 
gas can fall into the potential well before the radiation is converted, the more efficient the 
process, hence the appealing nature of compact objects. For black holes approximately 
Ae ~ GM/R^s ~ O.lc^ = 10^° erg g^^ (where R^g is the radius of the marginally 
stable orbit) can be released, and even more if the hole is endowed with a large angular 
momentum. This efficiency is over a hundred times that traditionally associated with 
thermonuclear reactions (Hydrogen burning releases 0.007c^ ~ 6 x 10^^ erg g~^). Since 
the 1960s, we have also learned about stellar-mass black holes, where M ~ 5 — 10 Mq, in 
Galactic binary systems and ultra-luminous X-ray sources which, with less confidence, 
we also associate with black holes, primarily on energetic grounds. 

In these objects, accretion (and the accompanying radiation) is usually thought 
to be limited by the self-regulatory balance between Newtonian gravity and radiation 
pressure. A fiducial luminosity is the Eddington limit associated with quasi-spherical 
accretion, at which radiation pressure balances gravity. If Thomson scattering provides 
the main opacity and the relevant material is fully ionized Hydrogen, then this luminosity 
is 



o"T being the Thomson cross-section. This may be converted to a mass accretion rate if 
one considers that the accretion luminosity is L-^dd = -^acc = GMM/R, giving 



where i?* is the radius of the compact object. Now clearly this applies strictly only in 
a quasi-spherical configuration and a steady state. Relaxing these assumptions allows 
for greater luminosities in transient events, or for configurations in which the energy 
release is somehow coUimated. However they will not be greater than the expression 
given above by more than a factor of a few. 




(1) 
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The photon luminosity, for the duration of a typical short burst (a few seconds at 
most), is thousands of times larger than that of any active galactic nucleus (thought 
to involve supermassive black holes), and is 12 orders of magnitude above the limit 
([1]). The total energy, however, is not very far off from that of other phenomena 
encountered in astrophysics, and is in fact reminiscent of that released in the core of a 
supernova. The Eddington photon limit ([1]) is circumvented if the main cooling agent 
is emission of neutrinos rather than electromagnetic waves. The associated interaction 
cross section is then many orders of magnitude smaller, and the allowed accretion rates 
and luminosities are correspondingly higher. For example, using the cross section for 
neutrino pair production, the Eddington limit can be rewritten as 

i^Edd,. = 8 X 10^3 {^^) ' (M/Mq) erg s'S (3) 

with an associated accretion rate, assuming unit efficiency for conversion of mass into 
neutrino energy, 

MEdd,. = OA{M/M^) [^^) ' Me s'^ (4) 

The time it would take an object to radiate away its entire rest-mass energy in this way 
is a mass-independent Eddington time given by 

M ^ . / E, ^2 



tKdd.u — —■ ~ 2.5 



Mem., VSOMeVy 
while the time scale over which an accretion-driven source would double its mass is 
~ {L / LEdd,u)~^ X (efficiency)"^ x tEdd,u- The dynamical time scales near black holes are 
modest multiples of Rg/c, where Rg is the characteristic size of the collapsed object 

Rg = GM/c^ ~ 1.5 X W^{M/Mq) cm, (6) 

and are therefore much shorter than tEdd,u- A fiducial Eddington density, characteristic 
near the horizon when the hole accretes at the Eddington rate, is 

It should be noted that the typical Thomson optical depth under these conditions is 

rT ~ 4(id,.^g ~ 10^^ (8) 

and so, as expected, photons are incapable of escaping and constitute part of the fluid. 
For completeness, we can also deflne an Eddington temperature, as the black body 
temperature if a luminosity LEdd,u emerges from a sphere of radius Rg, 

TEdd. = ( J^^^] ^ ~ 5 X lO^i (M/Mo)~i/^ ( — ^) K, (9) 
^"^"^'"^ {AnRlasBj v / o; V50MeV; ' ^ ^ 



or 



/ p \ -1/2 

kT^dd,. - 45 (M/Me)-V4 (^^^ j MeV, (10) 
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and an Eddington magnetic field strength 



BEdd,u = I 1 ~ 3 X 10 



1/2 . P X -1 



Finally, for comparison, we define Tth as the temperature the accreted material would 
reach if its gravitational potential energy were turned entirely into thermal energy. It 
is given by 

^ GMrUn 19 , , 

Some authors use the related concept of the virial temperature, T^ir = Tth/2, for a 
system in mechanical and thermal equilibrium. In general, the radiation temperature 
is expected to be < Tth. In deriving the above estimates we have assumed that the 
radiating material can be characterized by a single temperature. This may not apply, 
for example, when a hot corona deforms the neutrino spectrum away from that of a 
cooler thermal emitter |84j . 

Although we have considered here the specific case of neutrino pair creation, the 
estimates vary little when one considers for example coherent scattering of neutrinos 
by nuclei and/or free nucleons (except for the energy scaling). Similar overall fiducial 
numbers also hold for neutron stars, except that the simple mass scalings obtained here 
are lost. 

It is thus clear from the above estimates that when mass accretes onto a (stellar- 
mass) black hole or neutron star under these conditions, the densities and temperatures 
are so large (p ~ 10^^ g cm~^, T ~ 10^^ K) that: (i) photons are completely trapped; 
and (ii) neutrinos, being copiously emitted, are the main source of cooling since they 
can mostly escape. This regime, which requires correspondingly large accretion rates, 
is termed hypercritical accretion, and was considered for SN1987A shortly after the 
explosion [85l ES] . Such high accretion rates are never reached for black holes in XRBs 
or AGN, where the luminosity remains well below the photon Eddington rate ([1]). They 
can, however, be achieved in the process of forming neutron stars and stellar-mass black 
holes during the collapse of massive stellar cores. Note that at sufficiently high accretion 
rates, the density reaches the threshold for optical thickness even to neutrinos, which 
cannot then simply stream out (this occurs at p ~ lO^^g cm~^). 

2.2. Bestiary 

The current view is that SGRBs arise in a very small fraction (approximately ~ 10^^) 
of stars which undergo a catastrophic energy release event toward the end of their 
evolution. One conventional possibility is the coalescence of binary neutron stars 
[niliaiHTllHHlEHlEOlEIllgSllgnilM^ Double neutron star binaries, such as the famous 
PSR1913+16 [95], will eventually coalesce due to angular momentum and energy losses 
to gravitational radiation. The resulting system could be top-heavy and unable to 
survive as a single neutron star. However, a black hole would be unable to swallow the 
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large amount of angular momentum present. The expected outcome would then be a 
spinning hole, orbited by a torus of NS debris. 

Other types of progenitor have been suggested - e.g. a NS-BH merger [961 EH 198] , 
where the neutron star is tidally disrupted before being swallowed by the hole; the 
merger of a White Dwarf (WD) with a black hole [99]; the coalescence of binary WDs 
[m EH M [M]; or accretion induced collapse (AIC) of a NS [Ml [103], where the 
collapsing neutron star has too much angular momentum to collapse quietly into a 
black hole. In an alternative class of models, it is supposed that the compact objects 
are contained within a globular cluster, and that the binary system will evolve mainly 
through hardening of the binary through three-body interactions |104l I105j or physical 
star-star collisions |106j rather than by pure gravitational wave emission. Finally, 
violent reconfigurations of the magnetic field on magnetars [53] may offer another 
possibility. Table [2] provides a summary of the various rate estimates for some of 
these possible SGRB progenitors. Aside from the rate of SNe events, the rate of 
SGRBs and plausible progenitors in Table [2] are highly imprecise. The main sources 
of uncertainty are related to supernova kicks, the mass ratio distribution in binaries, 
mass limits for black hole formation, stellar radii and common envelope evolution 



[T071 [TM fT09l [TTO [ml [TT21 [TT31 [m fTT5l [TO im 



Table 2. Estimated progenitors of SGRBs and their plausible rates [66l fTOTl [TOT l fTTT I 
I116| I122j (in yr~^ Gpc~'^). The rates for these various occurrences are plagued by a 
variety of uncertainties in the supernova explosion mechanism, stellar evolution, and 
binary evolution, and should be taken as an order of magnitude estimate only. The 
oberved rates of Type (la + Ib/c) events given here provide a generous upper limit. 

Progenitor Rate (z = 0) 



The rate of gravitational mergers among BH-BH and BH-NS binaries depends on 
the orbital periods after spiral-in, and hence on the ill-understood details of the process. 
A binary with total mass lOMiM©, mass ratio q, and period (in days) in a nearly 
circular orbit will merge in 



Thus for a coalescence to occur within a Hubble time, spiral-in must have reduced the 
orbital period to < 1 day. Since a black hole can have a mass exceeding that of its 
companion helium core, spiral-in need not occur, as opposed to the case of a lower 



NS-NS 
BH-NS 
BH-WD 
NS AIC 

WD- WD 
SN Ib/c 
SN la 
SGRBs 



1-800 
0.1-1000 
0.01-100 
0.1 - 100 
3000 
60000 
150000 
10(47r/17) 
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mass object like a neutron star (except in very non conservative mass transfer). If the 
companion has a strong wind, spiral-in might never begin. If it does, envelope ejection 
might be less efficient than that by a neutron star, since black hole accretion could 
have a very low efficiency. If spiral-in does occur, however, the subsequent evolution 
may be simpler than in the double neutron star case. Helium stars less massive than 
3Mq expand their envelopes dramatically during core carbon burning |123] . requiring 
a second spiral-in to make a NS-NS binary which could merge in a Hubble time. The 
more massive He cores which could leave black hole remnants do not expand much, so 
a second spiral-in is not needed. 

Thus BH-NS binaries form at rates comparable to the NS-NS rate, but the fraction 
which merge depends on the miasma of mass transfer and spiral-in which determine the 
final period distribution. Gravitational waveforms of any merging systems will allow 
the masses of the merging bodies to be accurately determined |124l [T25j , shedding light 
on the underlying physics, as well as on the minimum black-hole mass determined by 
post-collapse infall. The much more certain NS-NS mergers will also allow the mass 
distribution of neutron stars to be determined. 

How might such a progenitor generate a relativistic outflow or a sudden release of 
electromagnetic energy? We now address various possible routes. 

2.3. Metabolic Pathways 

It has become increasingly apparent in the last few years that most plausible SGRB 
progenitors suggested so far (e.g. NS-NS or NS-BH, WD-BH mergers or physical 
collisions, and NS AIC) are expected to lead to a black hole plus debris torus system 
(Figure [7]). In this case there is no external agent feeding the accretion disk, and thus 
the event is over roughly on an accretion time scale (which would be on the order of 
one second). A possible exception includes the formation from accretion |127j or from 
a WD- WD [M] or NS-NS merger [1281 HSU [130] of a rapidly rotating (in some cases 
very massive) neutron star with an ultrahigh magnetic fleld (Figure [7]). If there is an 
ordered fleld B, and a characteristic angular velocity u, for a spinning compact source 
of radius i?*, then the magnetic dipole moment is ~ BR^. General arguments suggest 
|131[ I132j that the non thermal magnetic-dipole-like luminosity will be oc B"^ R^uj'^ / , 
and simple scaling from these familiar results of pulsar theory require flelds of order 
10^^ G to carry away the rotational or gravitational energy (which is ~ 10^^ erg) in a 
time scale of seconds |133l [I34j. 

One of the most perceptive theoretical discoveries that was made about black holes 
was that, when they spin, a fraction of their mass can be ascribed to rotational energy 
and is, in principle, extractable |135j . This is most convincingly demonstrated by 
observing that there exist orbits of test particles with negative total energy, (including 
their rest mass), within the ergosphere. If an infalling plasma cloud is attached to 
magnetic fleld lines anchored at a large distance (e.g., in an accreting torus), and the 
fleld drags the cloud backwards relative to the rotation of the hole placing it on an orbit 
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Figure 7. Schematic scenarios for plausible SGRB progenitors: neutron star mergers 
(b,d); neutron star (e) or white dwarf (f) disruption by a black hole; white dwarf 
mergers (a) ; accretion induced collapse of neutron star (g) ; and a recycled magnetar (c) . 
The dominant production channel for each scenario is depicted, where MS denotes the 
primary main sequence star. The (rough) relative in-spiral times due to gravitational 
radiation for compact mergers are shown (not to scale). (d,e,f,g) eventually lead to the 
formation of a black hole with a debris torus around it. Physical star-star collisions or 
three-body gravitational interactions rather than gravitational mergers could lead to 
similar conditions. A second way that the progenitor system may liberate its binding 
and rotational energy is to develop a large magnetic field and act as a magnetar (a,b,c). 
Panels adapted from [126j . 

with negative energy, the work that has been performed can be thought of as energy 
that has been effectively extracted from the spin of the black hole [136j . 

The binding energy of the orbiting debris, and the spin energy of the BH are thus 
the two main reservoirs for the case of a black hole central engine: up to 42% of the rest 
mass energy of the torus, and 29% of the rest-mass energy of the black hole itself can 
be extracted for a maximal black hole spin. SGRB activity can be powered by the black 
hole only as long as there is interaction with the surrounding gas, which will probably 
need to be centrifugally supported (although even for relatively low angular momentum 
a considerable amount of energy can still be extracted |137j ). 

The angular momentum is quite generally a crucial parameter, in many ways 
determining the geometry of the accretion flow. Even a little rotation can make a 
big difference, breaking the spherical symmetry and producing accretion disks instead 
of radial inflow, as envisaged originally by Bondi [138]. If the gas has no angular 
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momentum and the magnetic field is dynamically unimportant, there will be essentially 
radial inflow. Spherical accretion onto black holes is relatively inefficient despite the 
deep potential well, because the gas is compressed, but not shocked, and thus cannot 
easily convert gravitational to thermal energy. The flow pattern changes dramatically 
if the inflowing gas has a small amount of angular momentum. The quasi-spherical 
approximation breaks down when the gas reaches a radius -Rcirc ~ l^/GM, where I is 
the angular momentum per unit mass, and if injection occurs more or less isotropically 
at large radii, a familiar accretion disk will form (as occurs in X-ray binaries and white 
dwarf binary systems). The matter will instead dissipate its motion perpendicular to 
the plane of symmetry and form a differentially rotating disk, the rotational velocity at 
each point being approximately Keplerian, and then gradually spiral inwards as viscosity 
transports its angular momentum outwards. 

If the emission process is very efficient, the disk is dynamically cold and 
geometrically thin, in the sense that locally, kT -C GMrrip/R and the pressure scale 
height H{R) -C R. If gas passing through a thin disk reaches a radius within which 
the internal pressure builds up - either because it is unable to cool in an inffow time or 
because the radiation pressure force is competitive with gravity - the disk will become 
geometrically thick, with H R. In a thick disk or torus the pressure provides 
substantial support in the radial as well as the vertical direction, and the angular 
momentum distribution (now as a function of height as well as radius) may be far 
from Keplerian. Both types of conffgurations have been studied extensively and the 
agreement in some cases between observation and theory is extremely good, assuming 
a fundamentally empirical recipe for angular momentum transport |139l I140[ I14H I142j . 

As mentioned above in §2.11 in principle ffow onto a compact object can liberate 
gravitational potential energy at a rate approaching a few tenths of Mc^, where M is 
the mass inffow rate. Even for such high efficiencies the mass requirements of the more 
luminous SGRB sources are rather high, with 



where e is the overall efficiency. The inner regions of disks with mass ffuxes in this range 
are generally able to cool by neutrinos on time scales shorter than the inffow time. If 



a few gravitational radii in size, and the physical conditions can be scaled in terms of 
the Eddington quantities deffned in § 12.11 The remaining relevant parameter, related 
to the angular momentum, is finfiow/^^frecfaih where ffreefaii — {2GM / RY^"^ is the free fall 
velocity. The inward drift speed finflow would be of order ffreefaii for supersonic radial 
accretion. When angular momentum is important, this ratio depends on the mechanism 
for its transport through the disk, which is related to the effective shear viscosity. For 
a thin disk, the factor (fmfiow/^^freefaii) is of order a{H/Ry, where H is the scale height 
at radius R and a is the phenomenological viscosity parameter. 

Suppose that a given accretion rate yields a luminosity with an efficiency 0.1. 
Then the characteristic density, at a distance R from the hole, with account of the effects 




(14) 



m = M/M-Edd,u < 1, then the bulk of the neutrino radiation comes from a region only 
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of rotation, is 



p ~ m{R/Rg) ^/^(t^inflowMreefalOpEdd,, 



(15) 



and the maximum magnetic field, corresponding to equipartition with the bulk kinetic 
energy, would be 



Any neutrinos emerging directly from the central core would have energies of a few 
MeV. Note that, as mentioned above, kT^^^u is far below the virial temperature 
kT^ir ^ mpC^{R/Rg). 

The flow pattern when accretion occurs would be then determined by the value of 
the parameters L^/LEdd,u, which determine the importance of radiation pressure and 
gravity, and the ratio tcooi/^dynamicai, which fixes the temperature if a stationary flow 
pattern is set up. The preceding general discussion of neutrino-cooled accretion flows 
thus provides a basis for general models. So far we have discussed the power output 
that might be generated by the accretion process, but we have made no attempt to 
describe in detail the flow of the accreting gas. A hint that its dynamics may not 
be straightforward is provided by the existence of the Eddington limit for neutrinos, 
clearly illustrating that for the high accretion rates expected in SGRBs, forces other 
than gravity may be important. In addition, in many cases (and probably in most) the 
accreting matter possesses considerable angular momentum per unit mass which it has 
to lose somehow in order to be accreted at all. The reader is referred to ^ for fuller 
details. We now concentrate briefly on the origins of the jets which provide the power 
for SGRB sources. 

Two ingredients are necessary for the production of jets: first, there must be a 
source of material with sufficient free energy to escape the gravitational field of the 
compact object; second, there must be a way of imparting some directionality to the 
escaping fiow. Our eventual aim must be to understand the overall fiow pattern around 
a central compact object, involving accretion, rotation, and directional outfiow but we 
are still far from achieving this. Most current works who have discussed outfiow and 
collimation have simply invoked some central supply of energy and material. A self- 
consistent model incorporating outfiow and infiow must explain why some fraction of 
the matter can acquire a disproportionate share of energy (i.e., a high enthalpy). A 
brief summary of the various metabolic pathways is presented in Figure [HI 

The neutrino luminosity emitted when disk material accretes via viscous (or 
magnetic) torques on a time scale At ~ 1 s is roughly 



for a canonical radiation efficiency of 0.1. One fairly direct solution is to reconvert some 
of this energy via collisions outside the disk into electron-positron pairs or photons 
[92| I143j . If this occurs in a region of low baryon density (e.g., along the rotation 
axis, away from the equatorial plane of the disk) a relativistic pair-dominated wind 




(16) 
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Figure 8. Metabolic pathways for energy extraction. Panel a: Energy released as 
neutrinos is reconverted via collisions outside the dense core into e^ pairs or photons (i). 
The neutrinos that are emitted from the inner regions of the debris deposit part of their 
energy in the outer parts of the disk (ii), driving a strong baryonic outflow. This wind 
may be responsible for collimating the jet. Panel b: Strong magnetic fields anchored in 
the dense matter can convert the binding and/or spin energy into a Poynting outflow. 
A dynamo process of some kind is widely believed to be able to operate in accretion 
disks, and simple physical considerations suggest that fields generated in this way 
would have a canonical length-scale of the order of the disk thickness (i) . Open field 
lines can connect the disk outflow and may drive a hydromagnetic wind (ii). The above 
mechanism can tap the binding energy of the debris torus, while a rapidly rotating hole 
could contain an even larger energy reservoir, extractable in principle through MHD 
couphng to the exterior (iii;iv). Panel b adapted from |f44j . 



can be produced. An obvious requirement for this mechanism to be efficient is that the 
neutrinos escape (free streaming, or diffusing out if the density is high enough) in a time 
scale shorter than that for advection into the black hole. The efficiency for conversion 
into pairs (scaling with the square of the neutrino density) is too low if the neutrino 
production is too gradual, so this can become a delicate balancing act. Typical estimates 
suggest a lower bound of L uP ~ 10 "^L^/ when the entire surface area emits close to a 
single temperature black-body. The efficiency may be significantly larger if dissipation 
takes place in a corona-like environment [84] . 

One attractive energy extraction mechanism that could circumvent the above 
restriction in efficiency is a relativistic magneto hydrodynamic (MHD) wind |145[ 1134] . 
Such a wind carries both bulk kinetic energy and ordered Poynting flux, and it is 
possible that gamma-ray production occurs mainly at large distances from the source 
[531 11331 11341 11461 11481 1149j . A rapidly rotating neutron star (or accretion disk) releases 
energy via magnetic torques at a rate 

where P is the spin period, and B is the strength of the poloidal field at a radius R. 
The last stable orbit for a Schwarzschild hole lies at a coordinate distance R = 6Rg = 
9{M/Mq) km, to be compared with Rg = 3/2{M/Mq) km for an extremal Kerr hole. 
Thus the massive neutron disk surrounding a Schwarzschild black hole of approximately 
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2Mq should emit a spin-down luminosity comparable to that of a millisecond neutron 
star. A similar MHD outflow would result if angular momentum were extracted from 
a central Kerr hole via electromagnetic torques jl47 ]. The field required to produce 
Lem > 10^^ erg s~^ is enormous, and may be provided by a helical dynamo operating 
in hot, convective nuclear matter with a millisecond period [53]. A dipole field of the 
order of 10^^ G appears weak compared to the strongest field that can in principle be 
generated by differential rotation (~ 10^^[P/1 ms]^^ G), or by convection (~ 10^^ G), 
although how this may come about in detail is not resolved. Note, however, that it 
only takes a residual torus (or even a cold disk) of 10~^ MQto confine a field of 10^^ G. 
Orbiting debris with such large magnetic seed fields and turbulent fluid motions will 
give rise to a plethora of electromagnetic activity (Figure [8]). The topic lies beyond the 
scope of this paper and we refer the reader to the excellent review by Blandford |144j . 

A potential death-trap for such relativistic outflows is the amount of entrained 
baryonic mass from the surrounding medium. For instance, a Poynting flux of 10^^ erg 
could not accelerate an outflow to F > 100 if it had to drag more than ~ IO'^Mq of 
baryons with it. A related complication renders the production of relativistic jets even 
more challenging, because the high neutrino fluxes are capable of ablating baryonic 
material from the surface of the disk at a rate |150j 



Assuming that the external poloidal field strength is limited by the vigor of the 
convective motions, the spin-down luminosity scales with neutrino flux as L^ind ~ 
Lem (X (X oc Ll^^ , whcrc fcon is the convective velocity. The ablation rate given 
in equation f|T9l) then indicates that the limiting bulk Lorentz factor F^ of the wind 
decreases as L~^. Thus the burst luminosity emitted by a magnetized neutrino cooled 
disk may be self-limiting. Mass loss could, however, be suppressed if the relativistic 
wind were somehow coUimated into a jet. This suggests that centrifugally driven mass 
loss will be heaviest in the outer parts of the disk, and that a detectable burst may be 
emitted only within a relatively small solid angle centered on the rotation axis. 

2.4- Observational Tests 

In the preceding sections, we have endeavored to outline some of the basic physical 
processes that are believed to be of most relevance to interpreting SGRB sources, not 
so much because we believe that there is strong evidence in favor of them, but instead 
because it provides a framework in which to discuss the observations and to demonstrate 
that, as yet, SGRBs pose no threat to conventional physics. Of course, we are conscious 
that most observations tell us less about the primary source than the about secondary 
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Thus a rest mass flux limits the bulk Lorentz factor of the wind to 
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reprocessing of this power in the circumburst and galactic environment - from which 
we can learn about the central engine only by a chain of uncertain inferences. Unless 
SGRBs are eventually found to be accompanied by telltale emission features like the 
supernovae of long-duration GRBs, the only definitive understanding of the progenitors 
will come from possible associations to direct gravitational or neutrino signals. 

Spectral investigations will nonetheless be important as probes of the basic 
energetics and microphysical parameters of relativistic shocks [271 [32] . We are more 
sanguine about the possibihty of deriving accurate particle densities from afterglow 
observations. However, there is the possibility that occurrence in an unusually low 
density environment could cause the external shock to occur at much larger radii and 
over a much longer time scale than in usual afterglows [15H I152j . The X-ray afterglow 
intensity could be then below the threshold for triggering. This may be the case for 
SGRBs arising from compact binaries that are ejected from the host galaxy into an 
external environment that is much less dense that the ISM assumed for usual models. 
Another possibility for an unusually low density environment, made up only of very high 
energy but extremely low density electrons, is if the GRB goes off inside a pulsar cavity 
inflated by one of the neutron stars in the precursor binary |153[I154] . Such bubbles can 
be as large as fractions of a parsec or more, giving rise to a deceleration shock months 
after the SGRB with a consequently much lower brightness that could avoid triggering 
and detection. The difference between the low-density and high-density environment 
cases could be tested if future observations of afterglows reveal a correlation with the 
degree of galaxy clustering or with individual galaxies. 

The association of SGRBs with both star-forming galaxies and with ellipticals dom- 
inated by old stellar populations, to some researchers, suggested an analogy to type la 
supernovae, as it indicated a class of progenitors with a wide distribution of delay times 
between formation and explosion. Similarly, just as core-collapse supernovae are discov- 
ered almost exclusively in late-time star-forming galaxies, so too are long GRBs. Indeed, 
a detailed census of the types of host galaxies, burst locations and redshifts should help 
decide between the various alternatives suggested in § 12.21 This is because if the pro- 
genitor lifetime is long and the systemic kick is small, then the bursts should correspond 
spatially to the oldest populations in a given host galaxy. For early-type galaxies, the 
distribution would presumably follow the light of the host [39]. In contrast, a neutron 
star binary could take millions of years to spiral together, and could by then, especially 
if given a substantial kick velocity on formation, have moved many kiloparsecs from its 
point of origin. The burst offsets would then presumably be larger for smaller mass 
hosts. As new redshifts, offsets, and host galaxies of SGRBs are gathered, the theories 
of the progenitors will undoubtedly be honed. 
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3. Compact Object Mergers and Accretion Disk Assembly 

SGRB activity manifests itself on many different scales. The simplest hypothesis, 
however, is that the central prime mover or trigger is qualitatively similar in all events. 
The primary energy may be reprocessed in different ways, depending on the details of the 
circumburst (and galactic) environment. A brief summary of the various evolutionary 
pathways that may be involved in producing a SGRB is given in §2.21 The important 
message is that most current alternatives involve either an accreting stellar mass black 
hole (neutron star) or a precursor stage whose inevitable end point is such an object. 
Furthermore, as we discussed in §2.11 a black hole (neutron star) embedded in infalling 
matter offers a more efficient power source than any other conceivable progenitor; it 
is our firm prejudice that SGRBs are energized by mechanisms that involve either 
neutron stars or black holes. For this reason, in this section, we turn out attention 
to the formation of the trigger : the actual astrophysical system that is capable of, 
and presumably powers the observed burst. It is beyond the scope of this article to 
describe completely all evolutionary pathways, and we shall confine our attention to 
binary mergers and physical collisions. The study of compact binaries, particularly 
those involving pulsars, has evolved dramatically over the past thirty years, and carries 
a great deal of historical baggage. Consequently, most theoretical work has been directed 
towards describing these binary encounters. Our discussion will reflect this bias0. There 
are eight sections: § 13. II gives a summary of observed compact object masses; § 13.21 gives 
an overall picture of the merger scenario as we best understand it now; then in § 13.31 
we describe early considerations of these events; this is followed in § 13.41 by a summary 
of dynamical stability considerations in such systems; the merger itself is considered in 
§ 13.51 the effects of General Relativity and the fate of the remnant core are addressed in 
§ 13.61 and § 13.71 respectively: finally, the outcome of compact object collisions is studied 
in §[331 

3.1. The Mass of the Progenitor 

The mass of the compact object dictates a characteristic length and luminosity scale 
for SGRB activity, so we must consider what we know (and what we do not) from 
observations of systems containing compact objects in a more leisurely state of affairs. 

Measuring stellar masses, although a considerable observational challenge, is not 
nearly as difficult as measuring their radii (which is why the equation of state at 
nuclear densities remains one of the primary open questions in compact object and 
nuclear physics studies). When they occur in binaries, accurate mass determinations 
are possible (Figure [9]), and for cases where there are two compact objects in tight 
orbits, the constraints can be truly spectacular (PSRB1913+16 |156] . PSRB1534+12 
[T57] . and PSRJ0737-3039A,B [I58l [I59]). Neutron star masses are tightly clustered 

% The reader is referred to the exceUent review by Rosswog |155j for an alternative but complementary 
description of these source models. 
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Figure 9. Measured masses for neutron stars and black holes in binaries. The average 
values are Mns — 1.49Af0 and A/bh = 8.27Mq (where the scatter is much larger). 
The data are taken from [169j for neutron stars, and [170] for black holes. 



around the Chandrasekhar mass, at IAMq. Those that deviate farthest from this value 
towards high masses are found in Low Masss X-ray Binaries [160] . where the slow 
but steady mass transfer from a companion has produced top-heavy objects (the most 
dramatic example being J075H-1807 at 2.IM0 [161] ). At such an extreme, it is possible 
that conversion to strange (quark) matter |162l I163[ 11641 1165] may have occurred (see 
Jaikumar et al. for a review |166j ). Incredibly, it is possible that observations of 
superbursts (unusually long Type I X-ray bursts) may be able to discern if this is 
actually the case [1671 1168] . 
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Regarding black holes, the scatter is considerably greater, probably reflecting the 
fact that a threshold mass is necessary to create one in a core-collapse event, but no 
obvious hard limit exists at high masses, and a degree of variation is expected |171j . 
The constraints are also less stringent in most cases, due to a combination of effects, 
but nevertheless a boundary can be drawn at roughly 3Mq, separating the most massive 
neutron star from the least massive black hole. This is also close to the maximum mass 
allowed for a neutron star in hydrostatic equilibrium in General Relativity, based only 
on the condition that the equation of state must preserve causality |172] . One would 
thus expect that encounters between objects of this type would have mass ratios of 
about unity if two neutron stars are involved, and 0.2-0.3 if one of each interacts. 

Although white dwarfs are not usually considered as direct progenitors of GRBs, 
one should keep an open mind about their possible role, particularly in the light of 
unexpected discoveries made by Swift concerning late time variability and the inferences 
about the parent population of SGRBs. They may be responsible for the violent birth 
of neutron stars through accretion induced collapse or the merger of two dwarfs in a 
binary [T73l [IH [JTSl [IHl [IIHl [IIS], and it is possible that this remnant could 
power a GRB [Til EH I100[ llOlj . The mass of a white dwarf is somewhat lower that 
that of a neutron star, but they are larger, by a factor of at least / ~ 10^, leading to a 
natural (dynamical) time scale that is greater than for neutron stars and black holes by 
a factor /^/^ ~ 3 x lOl 

3.2. Waltzing Couples 

As a stellar binary revolves, gravitational waves carry away its energy and angular 
momentum. The loss rates depend on the system mass, semi-major axis, a, and 
eccentricity, and rapidly increase as orbital decay progresses. In the stage where 
a ^ R^, where R^ is the typical stellar radius, the evolution can be computed to high 
accuracy in the weak-field limit, using post-Newtonian expansions for point masses 



[T80t [T8T1 [T82| [T83l flM [T85| fTM [T871 [T88l [T89] . The characteristic decay time for 



which, when apphed to the binary pulsars PSR1913+16 [95] and PSRJ0737-3039 [T90] . 
gives to ~ 300 Myr and 85 Myr, respectively. These waveforms, calculated for the 
last few minutes of the binary system's life, will be required as templates for detection 
with interferometric gravitational wave detectors such as LIGO |191j and VIRGO [192j . 
For more common binaries with much longer periods, the decay rate will be negligible 
and the emission may be observed with the space-based mission LISA. For systems 
reaching small separations, the decay rate will increase catastrophically and a collision 
will ensue, followed by relaxation to a final, nearly steady state in which, in many cases, 
only a black hole will remain. This is the ring-down phase, in which the emission of 
gravitational waves and the configuration of the fossil black hole can be studied through 



circular orbits of period P is given by 




(21) 
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perturbative methods on a Kerr background (see |193j for a review). A composite 
gravitational radiation signal is shown in Figure [10] as an example. The decay at early 
times {t < 0) was computed using the weak-field limit approximation for point masses, 
while the merger waveform itself (t > 0) is the result of a three-dimensional numerical 
simulation for the disruption of a neutron star by a point mass. The initially slow 
decay accelerates (while the amplitude and frequency of the signal increase) and comes 
to an abrupt end when the separation is of the order of the stellar radius, giving a 
complex signal which depends on the details of the neutron star structure. In principle, 
observation of such radiation and its frequency power spectrum could help constrain the 
mass-radius relationship for neutron stars |194j . 

The middle ground, during which a violent redistribution of angular momentum and 
energy takes place, will produce a burst of gravitational, neutrino, and electromagnetic 
energy, perhaps in the form of a classical GRB [195] . There is no approximate, 
analytical solution to be found in this regime, and multidimensional numerical 
simulations are clearly required to address the behavior of the system. The process 
is highly dynamical, with thermodynamics and emission processes mattering little, and 
gravitational dynamics ruling the evolution (with the possible exception of systems 
containing white dwarfs). 

Even before high-resolution numerical simulations were carried out in essentially 
Newtonian gravity, it was assumed that the dynamical merger would indeed result in 
the formation of an accretion disk with enough mass and internal energy to account 
for the energy budget of a typical GRB, either through the tidal disruption of the 
neutron star in the former, or post-merger collapse of most of the central core in 
the latter. Compact binaries were thus seriously suggested as possible progenitors for 
GRBs at cosmological distances, with possible maximum power in excess of 10^'^ erg s~^ 

ini[l2l[88l[89l[90l[92l[93l[9l]. 

The characteristics of the problem, and the qualitatively different nature of the 
questions that can be of interest (or addressed) has dictated a range of approaches 
when searching for a solution. On the one hand, there is the question of generation of 
gravitational waves, clearly deserving of substantial effort in its own right, and tackled 
as such by several researchers. On the other, the generation of an electromagnetic signal 
(as a GRB, or another observable transient) requires the consideration of a different kind 
of input physics, mainly in the domain of thermodynamics, the equation of state under a 
wide array of conditions, nuclear physics, weak interactions and magnetohydrodynamics. 
Efforts in the field have generally proeeded along these largely orthogonal axes, with 
little overlap between them until quite recently, and it is only now that numerical 
simulations with any realistic microphysics are beginning to employ the full machinery 
of general relativity. It is thus only natural that results midway along this path have 
not always been in agreement, and indeed are not in some aspects even now. The recent 
and spectacular advances in the successfuU evolution of black hole binaries in numerical 
relativity give us today a glimpse of what may be possible to study in greater detail in 
the coming years (see the recent overview by J. Centrella [196] ). 
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Figure 10. Gravitational wave signals for the final stages of spiral-in of compact 
binaries consisting of a black hole and a neutron star with initial mass ratio q = 
Mns/A^bh — 0.3 and adiabatic index F = 3 (left column [239j ) and F ~ 5/3 (right 
column [238]). In this example, the observer is directly above the orbital plane of 
the binary along the rotation axis. The rise in amplitude and frequency shows the 
characteristic accelerated phase of decay, followed by a rapid transition to a ring-down 
after the violent tidal disruption of the star at t « 0. In the bottom panels, details of 
the merger and tidal disruption waveform are shown, from whose spectrum one could 
in principle extract valuable information regarding the equation of state at neutron star 
densities. Note that the neutron star with the more compressible equation of state is 
fully and rapidly disrupted, resulting in the total disappearance of the gravitational 
wave signal, whereas the one with a stiff equation of state leaves behind a low-mass 
orbiting remnant, which produces persistent emission at late times (see also Figure [TT]). 
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3.3. First Inroads 

To our knowledge, the merger of compact objects (i.e., neutron stars and/or black 
holes) was first considered by James Lattimer and David Schramm [1971 1198j . They 
initially studied the dynamics and disruption of a neutron star by a black hole in a 
bound system, and the discovery of the Hulse-Taylor pulsar in 1974 [95] soon made 
it clear that such encounters were an inevitable, if rare in terms of the lifetime of a 
typical galaxy, consequence of binary stellar evolution pathways. The computations 
were carried out in full General Relativity, with certain simplifications to make them 
tractable and their interest originally laid in the possibility of such events ejecting 
neutron-rich material to the interstellar medium, which might subsequently undergo 
r-process nucleosynthesis and contribute to the observed abundances of heavy elements 
[1991 1200] . GRBs had only recently been made known to the astrophysical community 
[201] . and the likelihood of them being of Galactic origin made Lattimer & Schramm 
consider such events to be unrelated. In their words: "Although the frequency of these 
events is too small to be important as far as currently observed anti-neutrino, 7-ray, 
or optical bursts are concerned, enough neutron star material may be ejected to be 
of nucleosynthetic importance." |198] . Over thirty years after their pioneering work, 
coalescing compact binaries stand today as one of the leading runners in models for 
SGRBs. It is noteworthy that only now are detailed computer simulations of such 
collisions beginning to fully consider the effects of General Relativity in the dynamical 
merger phase of tidal disruption (John A. Wheeler's "tube of toothpaste" mechanism 

[snaEna]). 

3.4. Dynamical Stability of Close Binaries 

The stiffness of the nuclear equation of state, or equivalently, the compressibility, is a 
parameter derived from microphysics that has a predominant role in the final outcome 
of a merger. It is most easily expressed by writing the pressure-density relation as 

P oc (22) 

where F is the adiabatic index. Figure [TT] shows the effect of its variation on the structure 
of a star by plotting the density profiles of static, spherically symmetric configurations 
in hydrostatic equilibrium. For F ^ 1, the result is a nearly-constant density (a stiff 
equation of state), while for F 4/3 (lower values will not produce a stable solution 
in Newtonian gravity) the star is highly centrally condensed (a soft equation of state). 
It is immediately obvious that when placed in the vicinity of a massive companion, the 
star made of more compressible material will feel tidal effects to a much lesser degree. 
The issue of how matter responds to pressure at extremely high densities is of course 
far from being solved. For example, while the stiffness of ordinary matter may be quite 
high at p ~ pnuc, with F ^ 2 — 3, the presence of exotic condensates at supra-nuclear 
densities could soften the equation of state considerably [204] . leading to potentially 
different behavior. 
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Figure 11. Normalized radial density profiles for spherically symmetric configurations 
of a given mass in hydrostatic equilibrium, for different compressibilities (parametrized 
through the adiabatic index F). The varying degree of central condensation makes the 
stiffer stars (with higher values of F) more susceptible to tidal effects when placed in 
a binary system. 

The equilibrium configurations for an incompressible fluid, and their stability, were 
considered by Chandrasekhar [205] . both for single and binary stars, in combinations 
that could represent double neutron stars or neutron stars in orbit around point masses. 
This work was extended later by Lai and collaborators [2061 ISQH [2081 ISQH MS\ to 
compressible configurations by the use of polytropic relations as equations of state. In 
a series of papers, they showed that purely Newtonian hydrodynamic effects can de- 
stabilize close binaries if the equation of state is sufficiently stiff (F > 2). This can 
occur before contact, when the separation is typically ~ 3 stellar radii. The effective 
binary potential becomes so steep that a dynamical (i.e., on an orbital time scale) plunge 
occurs, forcing the two stars to come together. This is independent of any consideration 
of General Relativity, and clearly shows just how powerful the departures from point- 
mass behavior can be when the combination of a finite stellar radius and compressibility 
are considered. For illustrative purposes, we show in Figure [I2] the equilibrium curves 
of total orbital angular momentum for a Newtonian neutron star-point mass binary for 
two compressibilities, as well as the point mass result (which is simply Jeq(-R) oc i?^/^). 
This has important implications for the initiation of mass transfer from one star to the 
other, the occurrence of tidal disruption and the final configuration of the system once 
it comes into near contact (see e.g.. Chapter 4 in [211j ). 
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Figure 12. Total dimensionless angular momentum as a function of separation (in 
units of the neutron star radius) for equilibrium configurations of BH-NS binaries with 
mass ratio q = 0.31 in a Newtonian approximation. The neutron star is modeled 
as a non-spinning tri-axial Roche-Riemann ellipsoid (see |207| ) of varying adiabatic 
index (curves for F = 3; 5/3 are shown). Tidal effects increase the steepness of the 
effective interaction potential and produce a minimum in Joq(i?), indicating the onset 
of a dynamical instability. The effect of the equation of state is clear, with the 
less compressible neutron star encountering instability, which will induce a dynamical 
plunge, at a greater separation than its more compressible counterpart. For reference, 
the equilibrium curve for point masses is also shown. As expected, it is a monotonically 
increasing function of radius and exhibits none of the turning points characteristic of 
the appearance of an unstable branch. 



The compressibility also impacts directly upon the mass radius relation of a star. 
For polytropes, we have 

R^(xM^*, (23) 

so the star will expand upon mass loss (or contract) for F < 2 (or F > 2). Thus it 
can further overflow its Roche lobe once mass transfer begins, in which case the process 
itself is unstable leads to complete disruption, or retreat from the Roche surface and, 
in principle, shut off mass transfer in the absence of external driving (which will not 
occur, because gravitational radiation emission will remove orbital angular momentum 
continuously). In either case, the details of the dynamics will depend on the global 
adiabatic index at high densities. 

3. 5. Investigating the Merger Phase 

The numerical studies of coalescing binaries began in the 1980s and 1990s, with 
computations of the gravitational wave emission and determinations of the stability and 
dynamics of the system at separations comparable with the stellar radius under various 
simplifying approximations [SEl EIHl EH EIS EIEl EIH EUl EIHl [2201 EHl [2221 [2231 
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MM EHZl [223 M M> 1223 [22SI [222] • The thermodynamical evolution of the fluid 
during merger, with the potential for electromagnetic energy release and GRBs was also 
considered by various groups [2301 ESU [2321 [2331 [23 [SMI |235l [EH]. These studies used 
Newtonian or post-Newtonian gravity with additional terms included in the equations 
of motion to mimic gravitational radiation reaction (which is the ultimate agent driving 
the binary evolution), and in some cases approximations to the equations of General 
Relativity. 

Figure [13] shows the dynamical evolution of a merging neutron star-black hole 
binary, during which the neutron star (modeled as a polytrope) is tidally disrupted and 
an accretion disk promptly forms |236l 12371 12381 12391 [98] . The general conclusion of this 
body of work (which includes a mock radiation reaction force to account for gravitational 
wave emission) is that when the orbital separation is of the order of a few stellar radii, the 
star becomes greatly deformed due to tidal effects, and the system becomes dynamically 
unstable. For stellar components with relatively stiff equations of state (F > 5/2) the 
instability is due to a strong steepening of the effective potential because of tidal effects 
(Figure [TT]) . In the opposite limit, the mass-radius relationship for the star is such that 
mass loss leads to an expansion and Roche lobe overflow, further accelerating the process 
and leading to a runaway in which disruption occurs. All of the features anticipated in 
the work of Lai and collaborators concerning the stability of the binary are apparent, and 
are seen as well in the numerical simulations of merging binaries performed by Rasio and 
collaborators |225ll226ll240j . The main additional result in terms of the dynamics is that 
configurations that are stable in principle up to Roche lobe overflow are de-stabilized 
by the mass transfer process itself, and fully merge as well within a few orbital periods. 
It is also quite clear that in the presence of gravitational radiation reaction (even in 
its crudest approximation) the system enters a dynamical infall at small separations 
that no amount of (even conservative) mass transfer can revert. Qualitatively, the 
result is largely independent of the mass ratio, the assumed initial condition in terms 
of neutron star spin (both tidally locked and irrotational configurations were explored 
in [2361 12371 12381 [239] ) and the assumed compressibility. Quantitatively, the details are 
different, and are primarily reflected in the final disk mass and the gravitational wave 
signal (see, e.g., [241] ). 

After a few initial orbital periods, a black hole of M ~ 3 — 5Mq, surrounded by a 
thick and hot debris disk (M^isk ~ 0.01 — O.IMq) approximately 4 x 10^ cm across is all 
that remains of the initial couple (Figure [Hj). 

For double NSs the situation is essentially the same in terms of the orbital dynamics 
prior to merger. Synchronization during spiral-in is impossible because the viscosity in 
neutron star matter is too small |242[ 1243] . and so they merge with spin frequencies 
close to zero, compared with the orbital frequency (such configurations are termed 
irrotational). Contact occurs at subsonic velocities and a shear layer then develops 
at the interface. Modeling this numerically is extremely difficult, because the layer is 
unstable at all wavelengths and vortices develop all along it [244] . possibly amplifying 
the magnetic field to extremely large values very quickly [130] . The mass ratio in double 
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Figure 13. The tidal disruption of a neutron star by a black hole in a binary: Each 
column (top to bottom in time) depicts the interaction (computed in three dimensions) 
for a given value of the adiabatic index T, going from stiff to soft (left to right). 
Logarithmic density contours of density (spaced every quarter decade, with the lowest 
one at \ogp[g cm^'^] — 11), are shown in the orbital plane, with additional shading for 
r = 2; 5/3. The distance scale is indicated in each panel, and a black disk of radius 2Rg 
represents the black hole. All three cases initially had identical mass ratios, q — 0.31 
pkl[239] . 
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Figure 14. During the tidal disruption of a neutron star by a black hole in a binary, an 
accretion disk is formed from some of the material stripped from the neutron star. A 
sequence of meridional slices through a three-dimensional simulation is shown, plotting 
logarithmic density contours for the density (contours are spaced every decade, with 
bold ones at \ogp[g cm~'^] =8; 12). In the first panel, the tidally deformed star is still 
visible to the left of the black hole. By 4 ms, the disk is growing around the black hole, 
and a slice of the tidal tail (see Figure [T^]) is visible to the left of the disk. By 17 ms 
the disk is practically fully formed, and is draining onto the black hole on a much 
longer time scale, determined by angular momentum transport effects. This particular 
calculation used F = 5/3, and the box shown covers 400 x 200 km |239j . 

NS binaries is quite close to unity, so a fairly symmetrical remnant would be a natural 
and expected outcome. However, in general even small departures from unity can have 
important consequences in the inner structure of the final object: the lighter star is 
disrupted and spread over the surface of its massive companion, which can remain largely 
undisturbed |245[ 12261 1228] . The final configuration now consists of a supra-massive 
neutron star (i.e., one with more mass than a cold, non-rotating configuration could 
support) surrounded again by a thick shock-heated envelope and hot torus similar to 
that formed in BH-NS mergers. The center of the remnant is rapidly, and in many cases, 
differentially rotating, which can have a profound impact on its subsequent evolution 
(see below, § 13. 7p . 

By definition, in either case the merger entails the formation of an object rotating 
at the break up limit. Material dynamically stripped from a star is thus violently 
ejected by tidal torques through the outer Lagrange point, removing energy and angular 
momentum and forming a large tail. For double neutron star binaries there are two such 
structures, and obviously only one in a BH-NS pair. The mass (typically up to a few 
tenths of a solar mass), energy and composition of these structures can be more easily 
studied by Lagrangian particle based methods, as they are several thousand kilometers 
across by the time the merger is over, and have thus moved off the numerical domain 
in grid based codes. By the inherent asymmetry in BH-NS binaries (in form and mass 
ratio), the black hole can receive a substantial kick (up to 10^ km s~^ [97J) during 
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the process, in addition to that arising from gravitational wave emission. Some of the 
fluid (as much as a few hundredths of a solar mass) in these flows is often gravitationally 
unbound, and could, as originally envisaged by Lattimer &; Schramm, undergo r-process 
nucleosynthesis |227l 1246] . The rest will eventually return to the vicinity of the compact 
object, with possible interesting consequences for SGRB late time emission (see §|5]). 

It is important to note that the debris disks are formed in only a few dynamical 
time scales (a few milliseconds), which is practically instantaneous considering their 
later, more leisurely evolution (which will be detailed below). 

The possibility that a fraction of the neutron star in a BH-NS binary would survive 
the initial mass transfer episode and produce cycles of accretion was explored at long 
time scales by Davies et al. [247j . This was motivated by the fact that numerical 
simulations employing relatively stiff equations of state at high binary mass ratios 
[971 1238[ 1248] showed that this might actually occur. However, this is most likely not 
the case, because: (i) the required equation of state was unrealistically stiff; and (ii) 
gravity was essentially computed in a Newtonian formalism. More recent calculations 
with pseudo-Newtonian potentials [2491 [27] and the use of General Relativity |250[ 1251] 
consistently fail to reproduce this behavior for a range of compressibilities in the equation 
of state. Instead, the neutron star is promptly and fully disrupted soon after the onset 
of mass transfer. 

3. 6. Effects of General Relativity 

Performing calculations with full General Relativity is clearly necessary, but is not 
an easy task. More importantly, it is not even clear how to estimate its effects with 
simple analytical considerations. Computations of the location of the innermost stable 
circular orbit (ISCO) in black hole systems indicate that tidal disruption may be avoided 
completely, with the star plunging directly beyond the horizon, essentially being accreted 
whole in a matter of a millisecond [252] . This would preclude the formation of a GRB 
lasting 10 to 100 times longer. But we stress that, just as for stability considerations 
in the Newtonian case, even if the location of the ISCO can be a useful guide in some 
cases, it cannot accurately describe the dynamical behavior of the system once mass 
transfer begins. In pseudo-Newtonian numerical simulations |249] and post-Newtonian 
orbital evolution estimates |253j . even when a near radial plunge is observed, the star is 
frequently distorted enough by tidal forces that long tidal tails and disk-like structures 
can form. The outcome is particularly sensitive to the mass ratio q = Mns/^bh and 
initial General Relativity calculations showed that the spin of the black hole is also 
very important, with rotating BHs favoring the creation of disks |254] . Dynamical 
calculations of BH-NS systems in pseudo-Newtonian potentials that mimic General 
Relativity effects typically show that for mass ratios q ^ 0.25 it is possible to form 
a disk, although of lower mass than previously thought, M^isk ~ 1O~^M0 [2491 [27] . 
A long one-armed tidal tail is formed as in the Newtonian calculations, with some 
material being dynamically ejected from the system. Recently, Faber et al. |250[ 1251] 
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have presented their first results for the dynamical merger phase of a BH-NS binary 
in General Relativity for low mass ratios {q ~ 0.1). While the the equation of state 
is simple (a polytropic relation with F = 2 was assumed), the final outcome is clearly 
dependent on the complex dynamics of mass transfer. Even though at one point in the 
evolution most of the stellar material lies within the analytically computed ISCO, tidal 
torques transfer enough angular momentum to a large fraction of the fluid, producing 
an accretion disk with Mdisk — O.IMq by the end of the simulation at t 70 ms. The 
black hole horizon is a point of no return, but the innermost stable circular orbit is most 
definitely not. 

Extending our own work, and further motivated by the first accurate localizations 
of SGRB afterglows (starting with GRB050509B, [6]), we computed the binary evolution 
of BH-NS systems in a pseudo-Newtonian potential in an attempt to better estimate 
the circumstances under which a disk may form, and its specific configuration if it does 
|237[ [27] . The most important parameters are the disk mass and size, as these set the 
global energy, density, and time scales for its evolution. The calculations employed the 
same formalism and code as before |239j . substituting the pseudo-Newtonian potential 
of Paczyhski & Wiita [255] . and using irrotational binaries as initial conditions (i.e., 
ones in which the neutron star spin is negligible). We explored variations in both the 
mass ratio (g ~ 0.1 — 0.3) and neutron star compressibility (T ~ 5/3 — 2), and found 
that full tidal disruption takes place as the star approaches the black hole, with a torus 
forming around it. It contains, however, substantially less mass than in the Newtonian 
case, because a larger fraction is directly accreted by the black hole. Additionally, a 
large portion of the remaining material is on highly eccentric orbits in the long tidal 
tails, and, while it will return to the vicinity of the black hole, it does not constitute 
immediately part of the accretion disk. Overall the torus mass is at least one order 
of magnitude lower than previously estimated, and somewhat colder, since it has not 
been shock-heated to the same degree by self-interaction of the accretion stream (see 
Figure [13]). Rosswog |249j has performed high resolution calculations of this type as 
well, but with the use of a realistic equation of state, and finds similar evolution and 
outcomes. 

In the case of merging neutron star pairs, general relativistic calculations [256[ 12571 
1258] have progressed even further, going beyond the adoption of a simple polytropic 
pressure-density relation [2591 I260[ 1261] . The outcome is predictably complicated, 
and depends sensitively on the total mass of the system and the initial mass ratio. 
Asymmetric binaries tend to produce more massive disks than those with identical 
components. As mentioned earlier for the case of Newtonian calculations, a small 
deviation from a mass ratio of unity is significant, since it can lead to the nearly 
complete tidal disruption of the lighter component, while the more massive star is 
largely unaffected. Shibata & Taniguchi [260] find that the lower the mass ratio, the 
more massive the resulting accretion disk, reaching Mdisk = O.O3M0 for q = 0.7. Both 
PSR1913+16 and PSRJ0737-3039 have q > 0.9 [1561 [M [IHHl [I69], so this mass ratio 
may be unrealistically low, but when q ~ 0.9 a considerable disk containing ^ lO^^M© 
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Table 3. Remnant disk masses in compact mergers for a variety of equations of state 
and various gravity methods: Newtonian (N); Paczynski & Wiita (PW); and General 
Relativity (GR). 
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still forms. This can easily account for the 10^^ — 10^° erg required for a typical SGRB, 
simply based on its gravitational binding energy. 

We give in Table [3] a summary of estimates for the disk properties for different 
progenitors, based on the calculations of various groups described above. The 
characteristic masses, maximum densities and temperatures are Mdisk/M© ~ 10^^ — 
10^^, p ^ 10^° — lO^^g cm~^ and T ~ 10^° — lO^^K respectively, and provide a starting 
point for the more detailed disk evolution calculations we will describe below in § HJ 

3. 7. The Fate of the Central Core 

In the case of merging neutron stars, the ultimate fate of the central object is still 
unresolved, and depends on the maximum mass that a hot, differentially rotating 
configuration can support. It has been known for some time that the maximum 
allowed mass can be increased by the effects of rotation (see, e.g., |265j ). and in 
particular, differential rotation [266]. In addition, the post-merger core is certainly 
not cold, as shock heating (at least for a mass ratio of unity) can raise the internal 
energy substantially. Shibata & Taniguchi |260] find that this threshold for collapse 
can be estimated as Mthres — l-35Mcoid, where Mcom is the corresponding value for 
a non-rotating and spherical, cold configuration. Based on the recent observation 
of PSRJ0751+1807 [M], M^om > 2.IM0, so a total mass greater than ~ 2.83M0 
is required for prompt collapse to a black hole. If Mcom < ^core < Afthres, various 
mechanisms (e.g., emission of gravitational waves, redistribution of angular momentum, 
magnetic fields) could act to dissipate and/or transport energy and angular momentum, 
possibly inducing collapse after a delay which could range from seconds to weeks (this 
has also been suggested as a two-step process which would be responsible for a SN/GRB 
association in the case of long GRBs [267j ). 



The Progenitors of Short Gamma-Ray Bursts 



39 



Differential rotation, besides providing additional support against collapse, could 
amplify seed magnetic fields in the core to large values |130i 11281 1129] and turn the it into 
a powerful magnetar. Spin-down torques, magnetic field winding and flaring alone could 
then in principle power a GRB, independently of the presence of a torus-like accretion 
structure |145l [97]. Note that this would also be one ingredient of the expected outcome 
of accretion-induced collapse of a WD, or a double WD merger. The surrounding 
environment would be much different, clearly, and its effect on the generation of a 
relativistic fireball would certainly be important (and possibly catastrophic |177j ). 
The magneto rotational instability (MRI) may also operate in such an environment, 
transferring angular momentum to the debris disk fast enough to keep it from being 
swallowed by the core when (and if) it collapses |268[ 1269] 1270] . An alternative possibility 
is that the core may deform into a bar-like structure, provided the equation of state is 
sufficiently stiff |260j . Gravitational torques in the presence of this asymmetry may then 
transfer angular momentum to the orbiting debris and increase the mass of the disk. In 
either of these last two cases one would then have a system similar to that occurring in 
a BH-NS merger. In general, the details remain highly uncertain, but powering a GRB 
is still a possible evolutionary pathway (see Figure 21 in [260]). 

3.8. Colliding Compact Objects 

The collision of unbound compact objects, rather than their merger in binaries, has not 
received much attention in the context of GRBs. Janka & Ruffert [231] computed the 
interaction of two identical neutron stars on parabolic orbits, and found a post-collision 
environment so contaminated with baryons that they concluded it was not viable as a 
GRB central engine. In addition, they estimated the event rate within galaxies, and 
found it much too low to be of interest. We have reconsidered this problem [106] . and 
find that if interactions in Globular Clusters (GC) are taken into account, the rates can 
have an important effect on the production of SGRBs, if the outcome of the collision is 
a favorable one. Exactly how frequent a collision takes place depends on the velocity 
dispersion and the stellar density in a given GC, as well as the number density of GCs 
in a given galaxy |106j . 

A similar scenario, the tidal disruption of stars by massive black holes (where the 
mass ratio is typically q = 10"'') was considered as a mechanism for feeding AGN, thus 
accounting for their luminosity [2711 1272] . Simple initial estimates and more detailed 
analytical calculations |273[ 1274] were subsequently confirmed in their fundamental 
aspects through numerical simulations by various groups |275l 12761 1277] . 

The relevant quantity fixing the strength of the interaction is 



where and Mbh are the stellar and black hole mass, Rp is the pericenter distance 
and R^, is the stellar radius. When rj ^ 1 the tidal field is sufficiently strong to disrupt 
the star (the energy for this is ultimately extracted from the orbital motion). For 




(24) 
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the scenario considered by Rees |274j . the star moves essentially in a fixed background 
metric since g ^ 1 (a condition fully exploited in the numerical calculations later carried 
out), and the typical pericenter distance is up to one hundred times the stellar radius. 
Thus the star does not directly impact its companion, and the interaction is entirely 
gravitational. It was found that approximately half the mass of the star is ejected from 
the system (at greater than escape velocity) and the other half remains bound. Using 
the Keplerian relation 



for the orbiting material, where de is the specific energy, and the fact that the differential 
mass distribution with energy, dM/de, is practically constant, the rate at which mass 
returns to the black hole is given by dM/dt oc 

Much can be learned from these calculations, and applied to the case of colliding 
compact objects of comparable mass. Much is different as well, primarily because the 
mass ratio is of order g ^ 1 — 0.1, rather than 10"''. This means that the pericenter 
distance is now comparable to the stellar radius. The initial encounter may thus resemble 
more a direct impact, with immediate exchange of mass, and a corresponding alteration 
to the mass ratio on a dynamical time scale. A second effect, as for mergers, is related 
to the compressibility of the material, and how the star responds in radius once mass 
transfer begins. 

We have now computed the outcome of parabolic encounters between black holes, 
neutron stars and white dwarfs, in various combinations, to determine the structure of 
the remnant and investigate its viability as a possible GRB central engine. The obvious 
difference between systems containing a WD and those that do not is that the time scales 
are typically longer in the former by a factor (_Rwd/-Rns)^''^- It is therefore natural to 
expect accretion and possible flares at hundreds, or even thousands of seconds, rather 
than tens of seconds after the main interaction. Whether these manifest themselves 
in 7-rays at all is another matter, and one that cannot be easily tackled numerically 
yet. The reason is that the WD radius is typically a thousand times larger than that 
of its companion, and hence the dynamical range that needs to be properly resolved is 
correspondingly increased. 

We flnd that collisions occur in two stages (see Figure [T5l) . During the flrst pe- 
riastron passage, an outward tidal stream is stripped from the low-mass member and 
ejected at high velocity. This act binds the surviving core, which rapidly returns for 
a second passage and is entirely disrupted, forming a second tidal tail, as well as an 
accretion disk around the massive companion. The duration of the evolution up to this 
point is comparable to the final disruption in a binary (about 10 ms). The end result 
is the formation of a system similar to that occurring in a binary merger, in that it 
contains a massive body (with the same uncertainties in delay and pathway to collapse 
to a black hole as before) surrounded by a debris disk of a few tenths to a few hun- 
dredths of a solar mass. We also considered the case of a white dwarf colliding with a 
massive neutron star (or low-mass black hole). We were unable at this point to follow 




de 



1 



(25) 
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HS-NS BH-NS W D-BH 




Figure 15. The collision of compact objects, and their remnants: The columns (top 
to bottom in time) show parabolic encounters between two neutron stars with masses 
1.75Mo and 1.4M0 (left), a neutron star and a black hole with masses 1.4M0 and 
4.5M0 (middle), and a white dwarf and a massive neutron star with masses 0.5Mq 
and 2.5M0 (right). All neutron stars are modeled with adiabatic indices T = 2, and the 
white dwarf with index T = 5/3, appropriate for a low-mass configuration. The scale 
is indicated in each panel, with logarithmic contours of density shown in the orbital 
plane (spaced every quarter decade, with the lowest one at logp[g cm~^] = 10; 10; 
for each column respectively). The two tidal tails and the accretion disk are clearly 
visible in the BH-NS collision. Note the large difference in spatial and temporal scales 
for the case involving the WD. The estimated delay for the WD core to return to the 
vicinity of the BH is 1000 s. All encounters have a parameter rj — 1, a.s defined in 
equation ([24)) . 
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the entire evolution of the system, because of the longer time scales involved, but a 
large tidal tail is also formed at first encounter, with the surviving core being bound 
as a result. Already, however, a disk is apparent, being a thousand times larger than 
in the NS/BH encounter. A detailed dynamical study of this type of collision will be 
presented elsewhere |106j . 

Despite a variety of possible initial configurations, our best guess for the ultimate 
prime mover thus consists of a similar astrophysical object: a stellar mass black hole 
surrounded by a hot debris torus; moreover the overall energetics of these various pro- 
genitors differ by at most a few orders of magnitude, the spread reflecting the differing 
spin energy in the hole and the different masses left behind in the orbiting debri^. On 
the hypothesis that the central engine involves hyperaccreting black holes, we would 
obviously expect the hole mass, the rate at which the gas is supplied to the hole, and 
the angular momentum of the hole to be essential parameters. In this case there is no 
external agent feeding the accretion disk (excluding tidal tail interactions, which we will 
consider in detail in § [5]), and thus the event is over roughly on an accretion time scale. 
Any attempt to explain their properties thus requires that we now consider its evolution 
and associated energy release. 



A possible exception includes the formation of a rapidly rotating (in some cases very massive) neutron 
star with an ultrahigh magnetic field. 
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4. Neutrino Cooled Accretion Flows and Short Gamma-Ray Bursts 

The way in which gas flows onto an accreting object depends largely on the conditions 
where it is injected. In this section we address the general dynamical evolution of 
neutrino cooled accretion flows under the conditions expected to occur in SGRB central 
engines. There are flve sections: § 14.11 is devoted to general considerations and a 
review of previous work addressing neutrino cooled flows; § 14.21 gives a summary of the 
microphysical and thermodynamical conditions in the fluid and considers the relevant 
equation of state and cooling processes in the disk; detailed calculations in a pseudo- 
General Relativistic potential are presented in § 14. 3t the generation of disk-driven winds 
and their accompanying outflows is discussed in § 14. 4t flnally, stability considerations 
are addressed in § 14.51 



4.1. General Gonsiderations 

Despite all of the complications inherent in the study of neutrino cooled flows (and 
accretion flows in general), two main ingredients, easier to qualify, are crucial for 
understanding the life and death of the accreting source. The flrst is simply the total 
mass available in the system at the outset. This flxes the global energy scales, both 
thermal and gravitational, and allows for an approximate determination of the available 
energy, with 

[sM-J [oI^^J [w^) ^^'-^''^ 

Note that we have estimated the thermal and gravitational energies to be of the same 
order. This is different than for classical thin disks, because the hypercritical accretion 
flows we are to consider are born in highly dynamical situations where there is an 
enormous store of internal energy (this is true for binary mergers as well as for collapsar- 
type scenarios). The fundamental reason is that in either case, internal energy was the 
primary source of support prior to the catastrophic gravitationally-induced collapse or 
tidal disruption. 

The second is the nature of the mechanism feeding the inner disk with mass and 
energy as a function of time. In the standard picture following core collapse, normally 
envisaged for long GRBs |278[ [85l I279j . the inflowing mass supplies an accretion rate 
M ~ 10~^Mq s~^ for a time tf ~ 10 s and is thereafter fed by injection of infalling 
matter at a rate that drops off as t~^^^. For the case of merging compact objects, it is 
usually assumed that no external agent feeds the previously assembled disk, and thus 
that the event is over roughly on an accretion time scale. This picture can be altered by 
the presence of tidal tails formed with material stripped from the stars during the initial 
merger phase and ejected on eccentric orbits. We will address this later. Regardless of 
the tidal tails, the physical conditions and instantaneous structure within the flow are 
similar in both scenarios, and conclusions drawn from one case may be used to infer the 
situation in the other. 
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The equivalent Eddington luminosity for neutrinos fixes a rough upper bound for 
the power output (considering coherent scattering off free nucleons as a source of opacity) 
of LEdd,u — 10^^ erg s ^. Combined with the available energy in the disk, one finds a 
characteristic lifetime at this power level of tEdd,u ~ ^burst — 10 ms. Two important 
effects alter this estimate: (i) the power is in fact much smaller than L-Edd,u, by at least 
two orders of magnitude, and (ii) accretion onto the black hole drains the disk of mass 
and energy on a viscous time scale, t^, which can be shorter if the efficiency of angular 
momentum transport is high. 

Speculation about the structure and evolution of a massive accretion disk 
surrounding a new-born black hole began as soon as they were suggested to be the 
prime movers of GRB sources [881 EOl |92] • Despite the large theoretical uncertainties, 
it was clear from the start that indeed dense, rapidly evolving structures with huge 
accretion rates (on the order of one solar mass per second) were to be expected, with 
natural time scales (dynamical and viscous) that could in principle account for GRBs. 
Whether or not the resulting signal turns into a burst or not is still a matter of debate, 
but as Stan Woosley put it |280j "If the signature of a SMq black hole accreting stellar 
masses of material in a minute is not a gamma-ray burst, what is it?". The physical 
regime at the expected densities and temperatures places these systems alongside core- 
collapse supernovae in the bestiary of astrophysics, and possibly makes them even more 
energetic. 

The detailed structure of neutrino cooled accretion flows (NDAFs) has been 
investigated by a number of groups over the past few years. Initially only steady- 
state, azimuthally symmetric and vertically integrated solutions were considered, with 
increasing level of detail mainly in what concerns the thermodynamics and the equation 
of state IMH [2821 |283l [281], but also in the effects of General Relativity [285]. These 
studies indicated that the instantaneous power output for plausible accretion rates, if 
maintained, could account for the overall energetics of GRBs. Neutrino cooling, when 
initiated at the proper temperature, would be an efficient mechanism for the removal of 
internal energy, and the disk would thin rather abruptly. The ID solutions of Popham 
et al. [281j were matched in the inner regions quite nicely with the initial numerical 
results for collapsar simulations [279] . Narayan et al. |282j considered among other 
things, how the radius of matter injection would affect the flow, and concluded that 
it needed to be rather small in order for neutrino cooling to be efficient. Otherwise, a 
large fraction of the accreting mass at large radii would, being unable to get rid of its 
internal energy, simply be blown away in a large-scale outflow. Kohri & Mineshige [283] 
then considered microphysical effects in greater detail, both in terms of the equation of 
state and the cooling processes. Finally, the effects of General Relativity on the disk 
structure have been studied by Chen & Beloborodov [285] for stationary flows around 
Schwarzschild black holes. 

More recently, time-dependent calculations in one |286j . two [2871 1288] -assuming 
azimuthal symmetry- and three dimensions |289l 1290] have tackled the question of 
time dependence and instabilities. The trade-off between relaxing the assumption of 
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0-symmetry or not is reflected in the time interval than can be reliably modeled at 
high resolution: tgim — 50 ms in 3D vs. tgim — 1 s in 2D. The former allows for a full 
exploration of azimuthal modes and instabilities, while the latter may provide reliable 
physical estimates on time scales comparable to the duration of SGRBs. Both 2D and 
3D simulations have relied on initial conditions taken from 3D binary merger calculations 
|239[ 1291] , and as such are a natural extension of these models. It has become clear from 
these studies that the neutrino energy release is rather rapid. As long as the accretion 
(or viscous) time scale tacc ~ ^disk/Af is longer than the cooling time tcooi ~ -E'int/-^!/, 
the fluid will be able to radiate most of its internal energy. The power is thus maintained 
at a fairly constant level for up to 100 ms, then drops rapidly. However, the density in 
the disk remains fairly high for essentially tacc? which can be as long as a few seconds. 
It is thus in principle able to anchor strong magnetic fields capable of driving MHD 
fiows. Another important result is that for high enough accretion rates (or equivalently, 
densities), the inner most regions of the disk can become opaque to neutrinos. This 
has important consequences for the cooling time scale, since the internal energy cannot 
escape immediately, but must diffuse out. It is also important for the composition, 
since a negative radial lepton gradient, akin to that occurring above proto-neutron 
stars following core-collapse |292j . can be established. Despite the stabilizing influence 
of differential rotation |293j . we have found that the disk can become convectively 
unstable, which can have important consequences for the generation of magnetic flelds 
and nucleosynthetic products in a possible outflow. 

General Relativity, as already pointed out, not only plays a crucial role in the 
merger dynamics of compact binaries, but may also be important in determining the 
evolution and energy output from the accretion disks thus formed. The disk is small 
enough that a substantial amount of material lies in the vicinity of the marginally 
stable orbit, where the potential departs signiflcantly from its Newtonian form. The 
black hole itself may be rapidly rotating, depending on the previous binary history, 
and thus produce inertial frame-dragging. Finally, the emitted neutrinos are subject 
to light bending effects, which can alter the efficiency for annihilation, and thus the 
corresponding spectrum. Dynamical simulations of post-merger disk evolution have 
been performed without |287l I288j and with |289l I290j account of these differences, 
through the use of a pseudo-potential of the form proposed by Paczyhski & Wiita ^255] 
and Artemova et al. [294j for rotating black holes. These effects have been taken into 
account in the computation of the energy deposition rates in the vicinity of the disk 
[93| mi I295[ I296[ I297[ 1298] , taking realistic configurations derived from compact mergers 
as input conditions in recent computations [299] . 

4-2. Physical Conditions and Relevant Processes 

These disks are typically compact, with the bulk of the mass residing within 4x 10'' cm of 
the black hole (which contains 3— 5M0). They are dense (10^ g cm"'^ ^ P ^ 10^^ g cm~^) 
and hot (10^ K < T < 10^^ K). The temperature is in fact high enough that the 
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nuclei are practically fully photodisintegrated in the inner regions. Neutronization 
due to the high densities can cause the electron fraction to drop substantially below 
Ye = 1/2, and capture of pairs by free neutrons and protons provides the bulk of 
the cooling for p > 10^° g cm~'^. The gas is composed essentially of a particles and free 
nucleons in nuclear statistical equilibrium (NSE), pairs (the relative importance of 
positrons is quite sensitive to the degeneracy parameter rj^ = Pe/kT, where /ig is the 
electron chemical potential), trapped photons and neutrinos of all species (e^ captures 
and annihilation produce only electron neutrinos, but nucleon-nucleon bremsstrahlung 
and plasmon decays can produce /i and r neutrinos as well). Studying the thermal 
evolution and the corresponding energy release thus requires a detailed equation of 
state and consideration of weak interaction rates and emission processes. In addition, 
the background gravitational field is intense enough that relativistic effects can come 
into play, at least in the inner regions of the flow. Finally, the whole situation can hardly 
be considered to be in a steady state, since it originates, in most cases, either from the 
collapse of the Fe-core in a massive star, or from the merger or collision of two compact 
objects. 

We have thus considered an equation of state in which the total pressure is given 

by 

P = Prad + Pgas + + P., (27) 

where is the pressure associated with neutrinos, Prad = aT'^/S, Pgas = (1 + 
3X^^c)pkT/ (4mp) and Xnuc is the mass fraction of photodisintegrated (and ideal) nuclei. 
As an approximate (but quite accurate) solution to the equations of nuclear statistical 
equilibrium between free nucleons and Helium (we do not consider the creation of iron- 
like nuclei) we use 

/ 2- x9/8/ \-3/4 / loj^x 

The degeneracy parameter of Fermions can vary over a large range in the disk, so it is 
necessary to use an expression (due to |300j ) that allows for it (but under the condition 
of relativity, which translates to p > 10^ g cm~^), namely 

1 



p = 

" 127T{hc) 



(29) 



and 



pYe 1 

n_ — 



3 



r]l + r]y{kTf (30) 



for the electron fraction. This formula reduces to the appropriate limits when the 
temperature is low (kT <ti rj^, implying P oc p^^^ for a cold Fermi gas) and when it is 
high {kT 3> T]e, giving P oc for relativistic pairs). Note that this automatically 
takes into account the presence of pairs in the limit of an ultra relativistic gas, making 
it unnecessary to alter the factor 1/3 appearing in the term for radiation pressure. 
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There is an additional effect which alters the equilibrium composition of the fluid, 
related to the optical depth. We approximate this by following Beloborodov [301] and 
computing the electron fraction as 

4. 0.437 (^?^) (31) 

in the optically thin and mildly degenerate case {Q = [rrin — mp]c^), and 

exp (32) 



Ye ^\ kT 

in the optically thick case (this follows from setting the neutrino chemical potential 
equal to zero in the reaction e + p ^ n and using Maxwell-Boltzmann statistics for the 
non-degenerate nucleons). We use an interpolated fit weighted with factors involving 
exp(— Tjy) to smoothly transition from one regime to the other. 

We have included neutrino energy losses from pair captures onto free nucleons by 
the use of tables |302j and pair annihilation from fitting functions [303J. Additionally 
plasmon decays |232j and nucleon-nucleon bremsstrahlung |304] losses are considered, 
although their contribution to the total luminosity is negligible. The resulting cooling 
rates are accurate over a wide range of temperature and density, and are vary consistently 
with the changing composition and conditions within the flow (particularly for the two 
most important contributions of capture and annihilation). 

We approximate the optical depth of the material to neutrinos by computing the 
cross-section for scattering off free nucleons and a particles heavy nucleons. These 
provide the dominant contribution, and are given by 

1 f \^ 

and 

a. = ao(^)^4sin^^H^]^ (34) 

where (Tq = 1.76 x 10^'^'' cm~^, 9w is the Weinberg angle and E^, is the energy of the 
neutrinos (roughly equal to the Fermi energy Ep of the mildly degenerate electrons). 
The optical depth is estimated as = H/ly, where H is the local disk scale height and 
is the neutrino mean free path. We find that in the inner regions, ~ 10^ in high 
mass disks. 

The cooling is then suppressed by a factor exp(— Tj,) to mimic the effects of diffusion, 
so the total neutrino luminosity is 

= j P"^(^ff + ^plasmon + (jpair + ^cap) exp(-r^)(im. (35) 

The transition from a transparent to an opaque fluid occurs at r,^ ^ 1, which corresponds 
roughly to p ^ lO^^g cm~^ (this happens largely in the high-mass disks, with their low- 
mass counterparts remaining mostly transparent). Photodisintegration cools the gas 
at a rate gphot = 6.8 x lQ^^{dX^-acl dt) erg s~^ cm~'^, which is included in the energy 
equation. Finally, angular momentum is transported through the disk by shear stresses. 
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Figure 16. Meridional structure and evolution of hyperaccreting flows around a 
Schwarszchild black hole, computed with the pseudo-Newtonian potential of Paczynski 
& Wiita |255| . The density, p (g cm^^) and cooling rate through neutrinos q 
(erg g~^ s~^) are plotted 100 and 200 ms after the start of the calculation, for disks 
with an initial mass (left panels) and O.OIM* (right panels), where Af» = 0.3Mq. 
The a viscosity coefficient is 10^^. Note the different bounds for the color scale between 
the two cases, in density as well as cooling. 



and we use the a prescription to vary the magnitude of the effective viscosity (all terms 
in the stress tensor are included in the momentum and energy equations |153j ). 

Clearly, although we have attempted to give an accurate thermodynamic description 
of the gas, several approximations remain. In the first place, we have assumed that all 
reactions considered reach equilibrium. This is quite accurate for photodisintegration, 
but not necessarily true for weak interactions. Second, our prescription for neutrino 
optical depth effects is obviously a far cry from the proper Boltzmann transport used 
in one-dimensional SN calculations (or even diffusion), but we believe it qualitatively 
captures the nature of the transparent-opaque transition, at approximately the correct 
density. Finally, we have not considered energy deposition from neutrinos back to the 
gas through absorption or incoherent scattering. This could be quite important in the 
outer regions of the disk (see § 12.31 and 14. 4p and relevant for the baryon loading of a 
possible GRB-producing outflow. 

4.3. Dynamical Evolution in a Pseudo-GR Potential 

Here we present results for our own 2D calculations in azimuthal symmetry for the 
evolution of hypercritical accretion flows around a black hole following compact object 
mergers. The initial conditions are essentially the same as those in [288], which were 
taken from 3D simulations of black hole-neutron star mergers, and are evolved in the 
pseudo-Newtonian potential of Paczynski & Wiita [255], which reproduces the existence 
of a marginally stable orbit at = 6GMbh/c^ for non-spinning black holes. The 
standard disk mass initially is M* = O.3M0. This is probably close to the highest mass 
that will be found in such disks, and we have also investigated conditions in lighter disks. 
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Figure 17. Color-coded optical depth in the p-T plane at 50 and 100 ms for disks 
with initial mass M, (left panels) and 0.1M» (right panels). The curved solid line 
marks the photodisintegration threshold between a particles and free nucleons (by 
mass), and the straight solid line marks the electron degeneracy threshold, given by 
kT = 7.7(p/10"g cm~3)^/3 MeV. As the disk drains into the black hole, it becomes 
more transparent, cooler and less degenerate. The a viscosity coefficient is 10^^. 



with Mdisk = 0.1,0.01M*. The microphysics included is the same as in our previous 
calculations (and described above in § I4.2p thus allowing for a clear identification of the 
effects of strong field gravity. 

The spatial structure of the disk is shown in Figure [16l where the cooling is also 
indicated. The overall structure is similar to that seen in the Newtonian simulations, but 



The Progenitors of Short Gamma-Ray Bursts 



50 



the maximum density is reduced, and only in the high-mass disks and at early times is 
the presence of an optically thick region (along the equator, and close to the black hole) 
evident. The thermodynamic structure is plotted in Figure fT7l where the regions in the 
density-temperature plane occupied by the disk are shown, along with the optical depth. 
The flow remains almost entirely photo dissociated in the inner regions, thus giving rise 
to the intense neutrino emission from pair capture onto free nucleons. It is also clear 
that the electron gas is neither ideal nor fully degenerate, and finite-degeneracy effects 
must be considered to compute its evolution properly. 

In our Newtonian calculations we found that neutronization was important in the 
inner regions of the flow. This is still the case for the relativistic calculations, although 
the effect is reduced, particularly for low disk masses. Figure [18] shows one snapshot 
of height-integrated radial profiles of the electron fraction in the disk for different 
disk masses. The most massive disk achieves higher densities and neutronization is 
substantially enhanced. The negative gradient in Ye{R) at small radii was also observed 
in Newtonian calculations, where it was a consequence of the transition to the optically 
thick regime. In the relativistic case it occurs even for low-mass disks, and is related to 
the appearance of a plunging region, where the radial velocity rapidly increases, once 
the fluid reaches the marginally stable orbit (see the bottom panel in Figure [T8l) . The 
net effect is also to produce an inversion of the radial lepton gradient, indicating that 
convection is possible (differential rotation has a stabilizing effect on this, and the full 
Solberg-Hoiland criterion must be considered, see below). 

The accretion rate and neutrino luminosity are plotted as functions of time in 
Figure [19] for a = 10~^ and disk masses covering two orders of magnitude. For 
comparison, the result for high disk mass in the Newtonian case is also shown. Clearly 
the interval during which a large luminosity is maintained is lower in the relativistic 
case, simply because the disk drains more rapidly onto the black hole. The accretion 
rates scale roughly linearly with the disk mass, whereas the power output in neutrinos 
shows a slightly steeper dependence (it also depends sensitively on the temperature), 
particularly evident at late times. 

One possibility for the driving of a relativistic outflow and the powering of a GRB 
is uV annihilation, and another is the often-quoted magnetic outflow |147j . If we assume 
a 1% efficiency for the first case at L^, = 10^^ erg s~^ (which scales with the square of 
the neutrino luminosity), and using the present set of calculations, we infer a scaling 
with disk mass for the total energy release as 



This has the same dependence as in the Newtonian case, but with a reduced intensity, 
by about a factor of five. For the magnetic scenario, the total energy release (assuming 
the magnetic field energy density is in equipartition with the internal energy density 




(36) 





(37) 
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Figure 18. Top: Height-integrated radial profiles of the electron fraction Yg for disks 
with varying mass in a pseudo-Newtonian relativistic potential. Bottom: Meridional 
velocity magnitude, (w^ + viy^^, for the same calculation. The marginally stable orbit 
is located at log(i?[cm]) w 6.6. In the outer disk, the gas consists mostly of a particles, 
and Ye = 1/2. As the fluid nears the black hole it rises slightly, then becomes smaller 
as neutronization takes place. The trend is then reversed and increases again. The 
minimum is a sensitive function of the disk mass (or equivalently, the density) . 



again in the same way as in the Newtonian case, but also reduced (this time by about 
a factor twenty). 

A comparison of these scahngs and the actual data for the prompt emission of four 
SGRBs detected in 2005 is shown in Figure [201 Clearly neutrinos are probably not 
the best way to produce a burst lasting more than a few tens of milliseconds, unless 
the disk mass is quite high. The energy release is somewhat insensitive to the burst 
duration in this case. On the other hand, magnetic energy extraction seems like a 
plausible mechanism, and would, at a given disk mass (if this is a standard quantity in 
any sense), reproduce the trend that longer bursts seem to have a greater fiuence. It is 
also apparent that a reduced disk mass is less of a problem for magnetic mechanisms 
than for neutrinos as regards the energetics. 

4.4- Disk Winds and Outflows 

The outer regions and surface of the accretion flow are subject to various effects that 
can unbind substantial quantities of mass and drive powerful winds, with important 
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Figure 19. The neutrino luminosity decreases as the disk drains into the black hole 
on a viscous time scale. The top (bottom) panel shows M (L^) as a function of 
time for calculations with three different initial disk masses, covering two orders of 
magnitude (M* = O.SM©) in the pseudo-Newtonian potential of Paczyhski & Wiita. 
For reference, the curve labeled [N] is the result of a calculation evolved in a Newtonian 
1/R potential. The viscosity parameter is set to a = 10~^ in all cases. 

implications for nucleosynthesis ^150[ \308\ 1309] . GRBs and supernovae |310j and the 
possible coUimation of relativistic outflows [143^ I31H I307[ 1312] . One of these is energy 
deposition by neutrino heating, another is thermonuclear burning and the corresponding 
energy release. The former has been studied extensively by a number of authors, in 
the supernova as well as GRB contexts. Dynamical disk evolution calculations do not 
always include this explicitly (our own simulations fall in this category), but reasonable 
estimates can nevertheless be derived in some cases. We include the latter in its crudest 
approximation, considering only the transition from a particles to free nucleons, whose 
mass fraction is given in equation (!28l) . 

Now the gravitational binding energy per nucleon is Eg^ ~ GM-Burrip/ R ~ 
5(10^cm/i?)MeV, and the production of one a particle releases 7.7 MeV/nucleon into 
the flow, so enough energy is available to produce a large-scale wind where the binding 
energy becomes sufficiently small. This is in fact seen in our calculations (Figure [2T|) . 
where the outward (mostly vertical) velocity is f„ ~ c/10 and Mwind — 2 — 5 x 
1O~^M0 s~^. The corresponding power is L^ind = ^wmd'i'wmd/2 — 1 — 3 x 10^° erg s~^ 
for a disk with M^isk — 0.3Mq initially. The mass outflow from a neutrino-driven wind 
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Figure 20. Top: Histogram of observed burst durations, taken from [305) . Bottom: 
comparison of energy vs. duration for GRB959509B (for which the dependence on 
redshift is indicated by the gray hne), GRB050709, GRB050724, and GRB051221A 
with estimates from compact binary mergers. The connected squares and circles 
show the total isotropic energy release (assuming coUimation into a solid angle 
51 — 47r/10) and duration (tgo) for vV annihilation and Blandford-Znajek powered 
events, respectively, as computed from our 2D disk evolution calculations in a pseudo- 
Newtonian potential |255j . The range in initial disk masses covers two orders of 
magnitude, and the effective viscosities are (from left to right), a — 10^^, 10^-^, 10^'^. 
The stars are estimates from z/j7-driven outflows in double neutron star mergers 
p6l[807] . 



in spherical symmetry was estimated by Qian & Woosley |150j as 

/ I \ 5/3 

Mwind ^ 5 X 10-^ " Mq s-\ (38) 

ylO^"^erg s J 

so at the disk neutrino luminosity of L^, ~ 10^^ erg s~^, it is interesting to note that the 
power of each mechanism is comparable (Table Hj). 



4-5. Stability and Convection 

An important question which cannot be addressed fully through steady-state 
calculations is that of stability. Several factors need to be considered, we address each in 
turn. The "runaway radial" instability discovered by Abramowicz et al. |313j applies to 
configurations with constant specific angular momentum as a function of radius. This is 
fully an effect of General Relativity, dependent on the existence of an inner Roche lobe 
for disk accretion, akin to the Li Lagrange point in binary stellar evolution. Essentially, 
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Figure 21. The energy released from He synthesis m the outer regions of the 
hypercritical accretion flow surrounding the black hole is sufficient to unbind the gas 
and produce a strong wind. The meridional velocity field shows an accelerated outflow 
in the regions where the mass fraction X of free nucleons drops appreciably (a threshold 
of 0.9 has been set in the plot), indicating the creation of a particles. 



Table 4. Energetics of winds. 
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if the fluid overflows this lobe, two things occur: (i) the black hole mass increases, 
pushing this critical point outward and (ii) material of the same angular momentum 
value is now in a position to be accreted. The effect is a runaway, in which the disk 
is accreted in a matter of a few dynamical time scales. It turn out that even a small 
positive gradient of the specific angular momentum, i{R), is enough to stabilize this 
condition, because matter with greater angular momentum will not be easily accreted, 
and this damps the instability quite rapidly. However, none of the post-merger accretion 
disks obtained through merger calculations have angular momentum distributions that 
are even remotely close to being constant, and so this turns out to be most likely 



The Progenitors of Short Gamma-Ray Bursts 



55 




t=200ms 



Figure 22. Convective stability of a neutrino-cooled accretion flow around a black 
hole with A/bh — 4Mq: The dark (light) areas represent regions where the Solberg- 
Hoiland criterion implies instability (stability) . The largest part of the unstable region 
clearly coincides with the optically thick region of the flow. As the disk drains into the 
black hole, the disk becomes more transparent and this region becomes smaller. The 
filament-like regions of instability close to the equator are due to the irregularities in 
the flow. These particular snapshots are taken from run a2M in [288], with a Shakura- 
Sunyaev viscosity coefficient a ~ 0.01 and 0.3 Mq in the disk initially. 

irrelevant for the evolution of GRB central engines. Likewise, the self-gravity of the 
disk is probably unimportant, particularly if one considers the results for disk masses 
from General Relativistic calculations to be more realistic. This has been raised as 
a possibility concerning the late-time X-ray flares in a few of the SGRBs discovered 
recently, and will be addressed in the following section. 

As for the standard photon-cooled thin disk solution, neutrino cooled flows were 
shown to be thermally unstable if radiation pressure dominates |282l I283j . Our numerical 
solutions which fully resolve the vertical as well as radial structure of the disk in time- 
dependent fashion show this to be the well. The difference is that the disks' 
response to this condition is a substantial thickening in the opaque regions, because of 
(i) the rise in pressure and (ii) the suppression of cooling. Rather than having a thin 
disk, the scale height H = \P / {dP / dz) \ ~ R. 
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We have also recently found that a convective instability is present in neutrino 
cooled accretion flows at high densities (or equivalently, accretion rates). This is 
analogous to what occurs above proto-neutron stars following core collapse |292] , and is 
directly related to the fact that at p ~ lO^^g cm~'^ the fluid becomes opaque to neutrinos 
(i.e., the equivalent optical depth, t^, is of order unity). A further rise in density is thus 
accompanied by an increase in entropy per baryon, Sf,, and in equilibrium electron 
fraction Y^. The classical Solberg-Hoiland requirement for convective stability can be 
written as two simultaneous conditions (e.g., [293] ): 



where P is the pressure, T is the temperature and s is the specific entropy. 

We find that the first of these is satisfied over most of the disk volume in dynamical 
calculations, and marginally so in the inner, opaque regions. Recall that established con- 
vection will erase the conditions which led to its occurrence in a characteristic turnover 
time icon — ^con/'?^con, uuless some external condition tries to maintain the instability, in 
which case (for efficient convection) marginal stability will ensue (see e.g., |314j ). The 
second condition is clearly not satisfied in the optically thick regime (see Figure [22]) and 
drives vigorous motions continuously, as can be seen from inspection of the correspond- 
ing velocity field (Figure 4 in [288] ). Although some circulations are also apparent in 
the optically thin region, these are of much smaller strength, and are essentially driven 
by the overshooting of fluid elements out of the convective region and into the outer 
disk, and their subsequent damping. 




(39) 



and 



1 dP Idf ds de ds 
p dz dR dz dz dR 



(40) 
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5. Prolonged Engine Activity 
5.1. Motivation 

It is possible that the mass and accretion rate shape many of the engine activity's 
visible manifestations. If we are to distinguish those properties of the activity which 
can be regulated by the engine itself and those which may varied independently, we must 
look for guidance from observations of the emitted radiation. A crucial observational 
development which must be noted is the discovery of late time X-ray flaring in a large 
number of bursts, both long and short [315[ 1316] . This has proven difficult to interpret 
in terms of refreshed shocks, because of their rapid rise and decay. A likely possibility 
is that it reflects renewed activity (although see |320j for an alternative view) in the 
central engine itself [3T71I3T81I3T9] . 

For short bursts, with durations of the order of half a second, these flares occur tens, 
or even hundreds of seconds after the main burst. There is also independent support 
that X-ray emission on these time scales is detected when light curves of many bursts are 
stacked [331 [M]. If the GRB itself is powered by a hypercritical accretion flow around a 
stellar mass black hole, the dynamical time scale is only milliseconds, while the viscous 
time scale, which is usually related to the burst duration, is at most a few seconds if the 
effective viscosity is equivalent to a ~ 10^^. It is thus in principle a problem to account 
for a resurgence of activity ten to one hundred viscous time scales (or equivalently, up 
to ten thousand dynamical time scales) later. 

Over the past year, a number of specific suggestions have been made concerning the 
production of fiares, both for long and short events - in some case, suggesting also that 
their occurrence in both types of events should indicate a common origin: interaction 
with a binary companion [1031 1321] ; fragmentation of a rapidly rotating core through 
non-axisymmetric instabilities [322] : magnetic regulation of the accretion fiow [323] ; 
fragmentation of the accretion disk through gravitational instabilities and subsequent 
accretion [324] : differential rotation in post-merger millisecond pulsars [38l 1325] and 
magnet ar-like activity. 

Despite all these suggestions, it remains unclear how exactly the regulation (in the 
case of magnetic mechanisms) or the fragmentation (in the case of self gravity) would 
come about. Magnetic regulation was actually addressed by van Putten & Ostriker 
|326] as a possible way to produce both long and short bursts from the same central 
engine, depending only on the spin of the black hole. Slow-spinning holes would be 
unable to halt accretion and thus a short burst would ensue. In rapidly spinning holes, 
the transfer of angular momentum through electromagnetic torques would stop or delay 
accretion, giving a long event. There are tantalizing clues in numerical simulations 
indicating that this might actually occur. Krolik et al. [327] find a large difference in 
the mass accretion rate onto the black hole in MHD simulations depending on the black 
hole spin. They report a calculation with a dimensionless Kerr parameter of a = 0.998, 
in which the accretion rate is strongly suppressed due to electromagnetic stresses, and 
no stationary state is achieved through the end of the computation (although they note 
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that the magnetic field geometry is fundamentally different from that suggested by van 
Putten & Ostriker). This effect is absent for lower values of the Kerr parameter. It is 
not clear how to extrapolate these results to the time scales involved in SGRB flares, 
however, since they occur at much later times (the calculations by Krolik et al. span 
about half a second for a black hole of four solar masses). 

5.2. Regulating the Accretion Flow 

One of the basic unknowns of accretion disk theory is the physical mechanism ultimately 
responsible for angular momentum transport and energy dissipation in the disk. It is 
well known that classical hydrodynamical viscosity cannot drive accretion at the rates 
inferred from observations in almost every astrophysical context where accretion disks 
are thought to play a crucial role. The usual way to overcome this difficulty is to assume 
that transport is dominated by some anomalous viscous phenomenon, possibly related 
to collective instabilities in the disk, and to give some adhoc parameterizations for its 
magnitude. 

It has been recently recognized that accretion disks threaded by a weak magnetic 
field are subject to magneto hydrodynamic instabilities (see Balbus & Hawley [328\ 
and references therein), which can induce turbulence in the disk and thereby transport 
angular momentum, promoting the accretion process. A possible alternative source 
of transport in cold disks is provided by gravitational instabilities |329j - although 
the outcome strongly depends on the thermodynamics of the disk. In particular, the 
fragmentation of a gravitationally unstable disk requires that the disk be able to cool 
very efficiently |330j . with a cooling time 



where is the local angular velocity in the disk. Gravitational instabilities would 
then lead to a self-regulating process: if the disk is initially cold, in the sense that 
kT ^ GM/ R, then gravitational instabilities would heat it up on the short dynamical 
time-scale, bringing it toward stability; on the other hand, if the disk is hot enough to 
begin with, radiative cooling will drive it toward an unstable configuration. As a result 
of these competing mechanisms, the switch associated with the onset of gravitational 
instabilities will act as a thermostat. For self-regulated disks, a simple relationship holds 
between the disk aspect ratio H/R and the mass ratio Mdisk/MeH- The disk is in fact 
marginally stable when its temperature is small enough that 



Of course, this relationship can only hold in the limit where Mjisk/MBH 1- 

Disk evolution calculations of neutrino cooled fiows have so far failed to show that 
the disk is close to instabihty [2831 [SHU [2881 [SHU [290], but this eventually may be 
the case when enough time has elapsed to allow the disk to significantly cool. It is 
not clear how to extrapolate these results to the time scales involved in SGRB fiares, 
since current calculations span at most a few seconds. It remains to be discussed 
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Figure 23. During tiie merger of a blaclt iiole and a neutron star, long one-armed 
tidal tails of material stripped from the star are formed. Their structure is plotted 
here for two values of the adiabatic index F of the neutron star fluid, with the scale 
indicated in each panel. The thick black line across the tail for F = 2 divides material 
that is bound to the central black hole from that which has enough energy to escape 
the system. Note the condensations at regular intervals in the tail for the stiff equation 
of state, produced by gravitational instability. 

under which conditions the development of a gravitational spiral structure would be 
indeed able to transport angular momentum in the disk efficiently, hence favoring 
accretion, and whether this might be able to account for flares. There is an additional 
concern with the the use of a viscous formalism in self-gravitating accretion disks: 
Balbus & Papaloizou [331] have shown that in general the energy transport provided by 
gravitational instabilities contains global terms, associated with wave energy transport, 
that cannot be directly associated with an effective viscosity. Their relative importance, 
and the conditions under which they may play a significant role are as yet poorly 
understood. 

5.3. Flares from Tidal Tails 

In this section we explore the possibility that the long tidal tails formed during compact 
object mergers and/or collisions may provide the prolonged mass inflow necessary for the 
production of late flares in SGRBs. The general idea of extended emission from compact 
object mergers has been advanced before [IT], relying on the gradual disruption of the 
neutron star core over a time scale that is much longer than the dynamical one |247j . 
As we have argued above, we do not believe this will actually occur, but find that the 
possibility of injecting matter into the central engine at late times remains, in a modified 
form. 

We noted already that large scale tidal tails are a common feature formed during 
mergers and collisions between compact objects. These are typically a few thousand 
kilometers in size by the end of the disruption event in the case of neutron stars and 
black holes, and a thousand times larger for events involving white dwarfs (see Figures [231 
and [24|) . The amount of mass contained in these structures is significant, and typically 
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Figure 24. During parabolic collisions of compact objects, similar structures as those 
found in mergers occur, (a) For a white dwarf-black hole collision, a large disk and 
tail form in a few minutes. The core of the white dwarf will return to the vicinity of 
the black hole in 10'^ s. (b) The collision of a neutron star with a black hole (the initial 
mass ratio is g = 0.31), using the pseudo-potential of Paczyhski & Wiita |255j . reveals 
the formation of a single large tail, with condensations due to self gravity appearing 
at late times (see Figure [23]). Note the different spatial and temporal scales between 
the two cases. 
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Figure 25. The differential mass distribution, dM/de with energy, computed for 
black hole neutron star mergers and one parabolic collision (labeled " C" ) . Each curve 
is marked with the initial mass ratio. The material with negative energies will fall back 
onto the central mass after following essentially ballistic trajectories. 

Mtaii ~ 0.01— O.O5M0. Depending on the details of the equation of state, a small fraction 
of this (10~^ — 1O~^M0) is actually unbound, and will escape to the surrounding medium. 
The rest will eventually fall back onto the leftover disk, and probably be accreted onto 
the central mass. 

Once the initial dynamical interaction is over, the fluid in the tails is practically on 
ballistic trajectories, moving in the potential dominated by the central mass. It is thus 
possible to compute the rate at which this matter will fall back and accrete onto the 
existing disk or black hole. To compare with the simple case of a supermassive black 
hole mentioned in § 13. 8[ we have also computed the differential distribution of mass 
with energy for this fluid: dM/de. Although it appears at first glance to be roughly 
constant, there are in fact significant deviations from this. The corresponding accretion 
rate thus differs from the t~^^^ law derived by Rees |274j . and is somewhat shallower, 
with Mfb oc The fall back accretion rate illustrated in Figure [26] indicates that 

the bulk of the matter in the tail returns to the black hole in roughly one second. 

Now, clearly one second is much too short to account for flares, which typically 
occur after 30-100 s. However, the fluid does not directly fall onto the black hole, since 
it has considerable angular momentum. In fact, its circularization radius is comparable 
to the size of the disk which was formed around the black hole in the first place (typically 
200-500 km). This is the point at which the gas will settle if enough of its energy is 
dissipated, circularizing its orbit. We thus have a situation in which in a short time 
scale, tfb ~ 1 s, a mass Mtao ~ O.OIM© returns to a radius i?circ ~ 5 x 10'' cm. If 
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Figure 26. The accretion rate onto tlie central region from the material in the tails is 
plotted for the same calculations as shown in Figure [25l Initially roughly constant, it 
rapidly assumes a power law decay, with index « 4/3. This is roughly independent of 
the type of interaction (although we note that at late times numerical noise makes an 
accurate determination of the slope more difhcult). The lower limit indicated at 100 
seconds marks the accretion rate required to account for the typical energy seen in a 
late X-ray flare in a SGRB. 



dissipation occurs (e.g., through shocks with what is left of the initial disk), a new ring 
of matter may form at R ^ Rare- The evolution of this newly injected mass will depend 
on the previous history of the original accretion disk, and on how much of it is left after 
the fall back time tfb. Evidently the component with the most energy will dominate the 
behavior of the system, and if Mdisk ^ Mt^n it is hard to see how the infalling tail could 
produce a substantial alteration of the overall flow properties and an accompanying 
observable signal. If the reverse is true and Mtaii ^ Mdisk, however, perhaps the inner 
accretion disk and the true accretion rate onto the black hole can be modified and a 
secondary episode of energy release be provoked. If we naively estimate the viscous time 
scale for the newly formed ring at i?circ as 

tvisc ~ RlJ lOu, (43) 

where 

u = acl/QKep (44) 

is the viscosity coefficient, Cg is the gas sound speed and ^Kep is the Keplerian angular 
frequency, then typical values from the calculations indicate that 
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Viscous time scales of 10 — 10'^ seconds are thus in principle plausible for 10~^ < 
a < 10~^, and are much longer than the injection interval itself. This is akin to the 
instantaneous injection of a discrete amount of matter at a given radius, and will produce 
an accretion episode with a delay proportional to its duration, tdeiay tacc (e.g., |211j ). 
naturally accounting for this observational fact in SGRB flares. 

An additional constraint comes from consideration of the total energy release in the 
observed flares. Taking neutrino emission as an example, the accretion efficiency is 

eacc = U/Mbuc' ~ 0.1, (46) 

based on our own calculations [2871 1288] and those reported by other groups [2891 1290j . 
where Mbh is the actual accretion rate that feeds the central black hole, and is reduced 
from the rate at which the disk itself is fed by a factor ( ^ 1/2, the remainder going 
into outflows (see also [282] ). The accretion efficiency is predominantly dependent upon 
the flow being optically thin to its own emission, and can thus be extrapolated down to 
the average accretion rate produced by the fall back material, 

Mfb ~ Mtaii/tvisc = 2 X 10-^Mq s-\ (47) 
Thus the neutrino luminosity would be 

Lu ^ eaccMBHc' = 3 X 10^^ erg s-\ (48) 
Accounting for a total fluence Lflare ~ 3 x lO^^erg s~^ and a total energy release 
Efinre ~ 3 X 10^^ erg over tfiare ~ 100 s (typical numbers for observed events) would 
require a uV annihilation efficiency of 10~^. This is close to being unrealistically high at 
this luminosity level, and would thus argue against neutrinos as an ultimate source for 
the flaring behavior. On the other hand, magnetic energy extraction could presumably 
also operate, since the injection can deposit enough mass in the disk to anchor sufficiently 
strong magnetic fields. It is hard to see, however, how either option could account for 
extremely delayed flares occurring after 10^ seconds, but perhaps those can be due to 
refreshed shocks [32] . 

5.4- Instabilities and Fragmentation in Fluid Tails 

The hydrodynamic stability of the ejected tail is a question that deserves special 
comment. The fluid moves globally and to a good approximation in the potential of 
the central object. It is, nevertheless, prone to an instability due to self-gravity, which 
may affect its structure on small scales. This is the varicose or sausage instability, and 
it can be shown (see Ch. XII in ref. |332j ) that for a cylinder of incompressible fluid, 
gravitational instability sets in for perturbations with wavelength 

A > A* = 2TxR^yx/x\ (49) 

where Rcy\ is the radius of the cylinder and x* ~ 1. The fastest growing mode, which 
determines the size of the fragments upon manifestation of this effect has x = 0.58, 
corresponding to a wavelength A = 10.8i?cyh and a growth time 

4 
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Binary merger calculations do not employ incompressible fluids, but for high enough 
adiabatic indices one can see that this behavior does indeed occur. 

The left panel in Figure [23] plots the large scale structure formed during a black 
hole-neutron star merger in a simulation with F = 3. The knots in the tail are spaced 
at regular intervals, with wavelength A ~ 2.5 x 10^ cm ^ 10i?taii, as predicted by the 
incompressible analysis. The estimated time scale for their formation is r ^ 20 ms, 
in good agreement with what is seen in the calculations. Binary neutron star mergers 
with stiff polytropic pressure-density relations performed by Rasio et al. [226] show ex- 
actly the same structures. Once fragmentation occurs, the individual clumps continue 
to move in the overall gravitational potential, and those that are bound to the central 
object will return to its vicinity at discrete intervals to inject matter. For lower values of 
the adiabatic index this behavior is no longer seen. Instead, the fluid expands smoothly 
over large volumes, and the tails are less well deflned. At the densities encountered 
in these structures, one would not expect a compressibility low enough to allow this 
instability to operate, and indeed it is not seen in merger simulations using realistic 
equations of state [249] . Nevertheless, it is worthy to consider their possible formation 
and subsequent evolution as a general feature. Note that if one substitutes the neutron 
star for a quark (strange) star, such droplets are a generic feature upon tidal disruption 
[333] . 
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6. Summary and Future Prospects 

In this final section we present a short summary of what we have learned so far about the 
physics of SGRB sources. Though the field is far from being mature, sufficient progress 
has been made in identifying some of the essential ingredients. We also describe the 
observational and theoretical prospects for the near future. 

6.1. Summary 

The progress in our understanding of SGRBs has been a story of consolidation and 
integration, and there is every indication that this progression will continue. Until 
recently, SGRBs were known predominantly as bursts of 7-rays, largely devoid of any 
observable traces at any other wavelengths. However, a striking development in the last 
several months has been the measurement and localization of fading X-ray signals from 
several SGRBs, making possible the optical and radio detection of afterglows, which in 
turn enabled the identification of host galaxies at cosmological distances. The presence 
in old stellar populations e.g., of an elliptical galaxy for GRB 050724, rules out a source 
uniquely associated with recent star formation. In addition, no bright supernova is 
observed to accompany SGRBs, in distinction from most nearby long-duration GRBs. 
It is now clear that short and long events are not drawn from the same parent stellar 
population. Even with a handful of SGRBs detected to date, it has become apparent 
that they are far from standard, both in their energetics and cosmological niche. This 
hints at the underlying possibility that the progenitor itself may be quite different from 
burst to burst, and not entirely restricted to the most discussed scenario involving the 
merger of compact binaries such as the Hulse- Taylor pulsar (although see |307j ) . 

The most fundamental problem posed by SGRB sources is how to generate over 
]^g50 gj,g ^Yie burst nucleus and channel it into collimated plasma jets. The cumulative 
evidence insistently suggests that the more powerful SGRBs must have "processed" 
upwards of IO^^Mq through a compact entity - only this hypothesis accounts for 
the high luminosity and compactness inferred from 7-ray variability. We believe that 
accretion onto a compact object, be it a neutron star or a stellar mass black hole, 
offer the best hope of understanding the "prime mover" in all types of SGRB sources 
although a possible attractive exemption includes a rapidly spinning neutron star with 
a powerful magnetic field. Consequently, most theoretical work has been directed 
towards describing the possible formation channels for these systems, and evaluate 
those which are likely to produce a viable central engine. The only way to cool the 
resulting hypercritical accretion fiow (other than through direct advective transport of 
the energy into a gravitational sink hole on a dynamical time scale) is by neutrino 
emission, circumventing the classical Eddington limit for photons and allowing for the 
conversion of gravitational binding energy into outfiowing radiation. Consideration of 
the associated energy release thus leads us down the path of detailed thermodynamical 
and microphysical processes unlike those encountered in most areas of astrophysics, 
except supernovae, where the physical conditions are quite similar in terms of density. 



The Progenitors of Short Gamma-Ray Bursts 



66 



entropy and internal energy. 

As of this writing, it is fairly clear from the concerted efforts of many groups working 
on different aspects of these problems that there is in principle no problem in accounting 
for the global energy budget of a typical SGRB from the class of systems here considered. 
The devil is in the details, of course, and the actual modes of energy extraction have 
yet to be worked out carefully. It would appear, however, that neutrino emission is 
more confined to be a competitive energy source in the early stages of the dynamical 
evolution, while magnetically powered events may be able to offer longer staying power. 

Still, various alternative ways of triggering the explosions responsible for SGRBs 
remain: NS-NS, NS-BH, BH-WD or WD- WD binary mergers, recycled magnetars, spun- 
down supra- massive NS and accretion induced collapse of a NS. Can we decide between 
the various alternatives? The progenitors of SGRBs are essentially masked by afterglow 
emission, largely featureless synchrotron light, which reveals little more than the basic 
energetics and micro physical parameters of relativistic shocks. In the absence of a 
supernova- like feature, the interaction of burst ejecta with a stellar binary companion 
|103j or with its emitted radiation [334j may be the only observable signature in the 
foreseeable future shedding light on the identity of the progenitor system. A definitive 
understanding will, however, come with the observations of concurrent gravitational 
radiation or neutrino signals arising from the dense, opaque central engine. 

Our understanding of SGRBs has come a long way since their discovery almost forty 
years ago, but these enigmatic sources continue to offer major puzzles and challenges. 
SGRBs provide us with an exciting opportunity to study new regimes of physics. 
As we have described, our rationalization of the principal physical considerations 
combines some generally accepted features with some more speculative and controversial 
ingredients. When confronted with observations, it seems to accommodate their gross 
features but fails to provide us with a fully predictive theory. What is more valuable, 
though considerably harder to achieve, is to refine models like the ones advocated here 
to the point of making quantitative predictions, and to assemble, assess and interpret 
observations so as to constrain and refute these theories. What we can hope of our 
present understanding is that it will assist us in this endeavour. 

6.2. Observational Prospects 

High energy astrophysics is a young field. It owns much of the remaining unexplored 
"discovery space" in contemporary astronomy. Two examples of this discovery space are 
at the extremes of observation of the electromagnetic spectrum. At the high end, there 
are already a few tens of TeV sources, while at the low end of < 50 MHz radio astronomy 
there are essentially no sources. Neutrino astronomy claims only two cosmic sources so 
far, the sun and SN1987a. Finally, as many of the most interesting high energy sources 
are ultimately black holes and neutron stars, the exciting field of gravitational wave 
astronomy — perhaps a decade away from birth — is inextricably linked to high energy 
astrophysics. These are the next frontiers and we consider the role of SGRBs within 
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them in turn. 

Atmospheric Cerenkov techniques are being used to detect 7-rays in the GeV-TeV 
range. These are important as both sources and as probes. Persistent sources have been 
identified with pulsars, blazars and supernova remnants and in each case are hkely to 
provide the best approach we have to understanding the fundamental nature of these 
sources. A tentative > 0.1 TeV detection of a long GRB has been reported with the 
water Cherenkov detector Milagrito [335] . Here the big question is to determine whether 
the SGRB jets comprise ultrarelativistic protons, that interact with either the radiation 
field or the background plasma, or if they are e^ pairs. TeV sources should be far more 
plentiful in the latter case. The combination of GLAST and telescopes like HESS and 
VERITAS ought to be able to sort this out. 

The same shocks which are thought to accelerate the electrons responsible for the 
non-thermal 7-rays in SGRBs should also accelerate protons |336l I338[ 1337] . Both the 
internal and the external reverse shocks are mildly relativistic, and are expected to lead 
to relativistic protons [339] . The maximum proton energies achievable in SGRB shocks 
are Ep ~ lO^*' eV, comparable to the highest energies measured with large cosmic 
ray ground arrays [340] . For this, the acceleration time must be shorter than both 
the radiation or adiabatic loss time and the escape time from the acceleration region 
[341] . The accelerated protons can interact with the fireball photons, leading to charged 
pions, muons and neutrinos. For internal shocks producing observed 1 MeV photons 
this implies > 10^^ eV protons, and neutrinos with ~ 5% of that energy, e,^ > 10^^ eV 
[342] . Another copious source of target photons in the UV is the afterglow reverse 
shock, for which the resonance condition requires higher energy protons leading to 
neutrinos of 10^'' — 10^^ eV |343j . Whereas photon-pion interactions lead to higher 
energy neutrinos and provide a direct probe of the shock proton acceleration as well as 
of the photon density, inelastic proton-neutron collisions may occur even in the absence 
of shocks, leading to charged pions and neutrinos [344] with lower energies than those 
from photon-pion interactions. The typical neutrino energies are in the ~ 1-10 GeV 
range, which could be detectable in coincidence with observed SGRBs. This is the 
province of projects like AMANDA, IceCube and ANTARES. Success in the former will 
suggest that ultrarelativistic outfiows comprise mainly protons. Neutrino astronomy 
has the advantage that we can see the universe up to ~ EeV energies. By contrast, the 
universe becomes opaque to 7-rays above ~ TeV energies through absorption by the 
infrared background. 

Finally the last and most challenging frontier is that of gravitational radiation, 
which is largely unknown territory. A time-integrated luminosity of the order of a 
solar rest mass (~ 10^^ erg) is predicted from merging NS-NS and NS-BH models, 
while that from collapsar models is less certain, but estimated to be lower. Ground- 
based facilities, like LIGO, TAMA and VIRGO, will be seeking such stellar sources. 
The observation the associated gravitational waves would be facilitated if the mergers 
involve observed SGRB sources; and conversely, it may be possible to strengthen the 
case for (or against) NS-NS or NS-BH progenitors of SGRBs if gravitational waves 
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were detected (or not) in coincidence with some bursts. The technical challenge of 
achieving the sensitivities necessary to measure waves from assured sources should 
not be understated. It may well take more than another decade to reach them. 
Neither, however, should the potential rewards. Gravitational waves offer the possibility 
of observing in an entirely new spectrum, not merely another window in the same 
(electromagnetic) spectrum. Furthermore, they will by their very nature tell us about 
events where large quantitites of mass move in such small regions that they are utterly 
opaque and forever hidden from direct electromagnetic probing, and are the only way 
(except perhaps for neutrinos) through which we can learn about them. There have been 
few regions of the electromagnetic spectrum where observation conformed to previous 
expectation. So it would be indeed remarkable if gravitational wave astronomy, or any 
of the other frontiers, turned out to be as we have described. Indeed, if past experience 
is any guide, they will surely provide us with new surprises, challenges and potential 
loose threads through which we can unravel another piece of the fabric of the universe. 

6.3. Theoretical Prospects 

Although some of the features now observed in GRB sources (especially afterglows) were 
anticipated by theoretical discussions, the recent burst of observational discovery has 
left theory lagging behind. There are, however, some topics on which we do believe that 
there will be steady work of direct relevance to interpreting observations. 

One of the most important is the development and use of hydrodynamical codes 
for numerical simulation of SGRB sources with detailed physics input. Existing two 
and three dimensional codes have already uncovered some gas-dynamical properties 
of relativistic flows unanticipated by analytical models (e.g., [345] ) . but there are 
some key questions that they cannot yet address. In particular, higher resolution is 
needed because even a tiny mass fraction of baryons loading down the outflow severely 
limits the maximum attainable Lorentz factor. We must wait for useful and affordable 
three dimensional simulations before we can understand the nonlinear development of 
instabilities. Well-resolved three dimensional simulations are becoming increasingly 
common and they rarely fail to surprise us. The symmetry-breaking involved in 
transitioning from two to three dimensions is crucial and can lead to qualitatively new 
phenomena. A particularly important aspect of this would be to link in a self-consistent 
manner the flow within the accretion disk to that in the outflowing gas, allowing for 
feedback between the two components. The key to using simulations productively is to 
isolate questions that can realistically be addressed and where we do not know what 
the outcome will be, and then to analyse the simulations so that we can learn what is 
the correct way to think about the problem and to describe it in terms of elementary 
principles. Simulations in which the input physics is so circumscribed that they merely 
illustrate existing prejudice are of limited value! 

A second subject ready for a more sophisticated treatment is related to the intensity 
and shape of the intrinsic spectrum of the emitted radiation. Few would dispute the 
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statement that the photons which bring us all our information about the nature of GRBs 
are the result of particle acceleration in relativistic shocks (e.g., [346] ) . Since charged 
particles radiate only when accelerated, one must attempt to deduce from the spectrum 
how the particles are being accelerated, why they are being accelerated, and to identify 
the macroscopic source driving the microphysical acceleration process. 

Collisionless shocks are among the main agents for accelerating ions as well 
as electrons to high energies whenever sufficient time is available (e.g., [3471 I348j ). 
Particles reflected from the shock and from scattering centres behind it in the 
turbulent compressed region have a good chance of experiencing multiple scattering 
and acceleration by first-order Fermi acceleration when coming back across the shock 
into the turbulent upstream region. Second-order or stochastic Fermi acceleration in 
the broadband turbulence downstream of collisionless shocks will also contribute to 
acceleration. In addition, ions may be trapped at perpendicular shocks. The trapping 
is a consequence of the shock and the Lorentz force exerted on the particle by the 
magnetic and electric fields in the upstream region. With each refiection at the shock 
the particles gyrate parallel to the motional electric field, picking up energy and surfing 
along the shock surface. All these mechanisms are still under investigation, but there is 
evidence that shocks play a most important role in the acceleration of cosmic rays and 
other particles to very high energies. 

There is no in situ information available from astrophysical plasmas. So one is 
forced to refer to indirect methods and analogies with accessible plasmas, found only 
in near-Earth space. Actually, most of the ideas about and models of the behaviour of 
astrophysical plasmas have been borrowed from space physics and have been refitted to 
astrophysical scales. However, the large spatial and long temporal scales in astrophysics 
and astrophysical observations do not allow for the resolution of the collisionless state 
of the plasmas. For instance, in the solar wind the coUisional mean free path is of the 
order of a few AU. Looked at from the outside, the heliosphere, the region which is 
affected by the solar wind, will thus be considered collision dominated over time scales 
longer than a typical propagation time from the Sun to Jupiter. On any smaller and 
shorter scales this is wrong, because collisionless processes govern the solar wind here. 
Similar arguments apply to stellar winds, molecular clouds, pulsar magnetospheres and 
the hot gas in clusters of galaxies. One should thus be aware of the mere fact of a 
lack of small-scale observations. Collisionless processes generate anomalous transport 
coefficients. This helps in deriving a more macroscopic description. However, small- 
scale genuinely collisionless processes are thereby hidden. This implies that it will be 
difficult, if not impossible, to infer anything about the real structure, for instance, of 
collisionless astrophysical shock waves. The reader is refer to [349j . and [350j for an 
excellent presentation of the basic kinetic collisionless (space) plasma theory. 

The most interesting problem remains, however, in the nature of the central engine 
and the means of extracting power in a useful coUimated form. In all observed cases of 
relativistic jets, the central object is compact, either a neutron star or black hole, and is 
accreting matter and angular momentum. In addition, in most systems there is direct or 
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indirect evidence that magnetic fields are present - detected in the synchrotron radiation 
in galactic and extragalactic radio sources or inferred in collapsing supernova cores from 
the association of remnants with radio pulsars. This combination of magnetic field and 
rotation may be very relevant to the production of relativistic jets (e.g., [iB]). Much 
of what we have summarized in this respect is conjecture and revolves largely around 
different prejudices as to how three dimensional flows behave in strong gravitational 
fields. There are serious issues of theory that need to be settled independently of the 
guidance we obtain from observations of astrophysical black holes. 
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